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Abstract: We study small charged black holes in global AdS spacetime in the presence 
of a charged massless minimally coupled scalar field. In a certain parameter range these 
black holes suffer from well known superradiant instabilities. We demonstrate that the 
end point of the resultant tachyon condensation process is a hairy black hole which we 
construct analytically in a perturbative expansion in the black hole radius. At leading 
order our solution is a small undeformed RNAdS black hole immersed into a charged 
scalar condensate that fills the AdS 'box'. These hairy black hole solutions appear in a two 
parameter family labelled by their mass and charge. Their mass is bounded from below 
by a function of their charge; at the lower bound a hairy black hole reduces to a regular 
horizon free soliton which can also be thought of as a nonlinear Bose condensate. We 
compute the microcanonical phase diagram of our system at small mass, and demonstrate 
that it exhibits a second order 'phase transition' between the RNAdS black hole and the 
hairy black hole phases. 
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1. Introduction 

There has been much recent interest in the physics of charged black brane solutions of the 
Lagrangian 1 
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1 We have chosen to work in 5 spacetime dimensions, and chosen the scalar field below to be massless 
merely for simplicity. The analysis of this paper carries over, without qualitative modification, in arbitrary 
spacetime and for arbitrary scalar potential. 
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where G5 is the Newtons's constant and the radius of AdSs is set to unity 2 . This system 
(sometimes called the massless Abelian Higgs model) admits a well known set of charged 
black brane solutions which are asymptotically Poincare AdS. Recent interest in this system 
is due to Gubser's observation Q that, at large e and when they are near enough to 
extremality, these black branes are unstable . The end point of the tachyon condensation 
sparked by this instability is a so called hairy black brane - a solution with a planar 
horizon immersed in a charged scalar condensate. Black branes interacting with such 
matter condensates are novel and interesting, and have been studied intensively over the last 
few years (see [§, ||], f|, ||, § and references therein) . Unfortunately almost all constructions 
of these solutions have been numerical 3 . 

Under the AdS/CFT correspondence, these planar AdSs solutions are dual to the 
states of a conformal field theory living on the flat spacetime R 3,1 . Another natural arena 
to study 3 + 1 dimensional conformal field theories is to work on 5 3 x i?timc instead. 
States of such a boundary field theory living on S 3 are dual to gravitational solutions that 
asymptote to global AdSs instead of planar AdSs. The corresponding charged black holes 
in global AdSs spacetime are characterised by their radius in units of the AdSs radius 
and their charge. At large horizon radius, these black holes are locally well approximated 
by black branes and we expect their physics to be qualitatively similar to the Poincare 
AdS charged branes. It is natural to enquire about the opposite limit: do small hairy 
black holes exist, and what are their qualitative properties? In this paper we answer this 
question by explicitly constructing a set of spherically symmetric hairy charged black holes 
whose radii are small compared to the AdSs radius 4 . Our construction is perturbative in 
the radii of our solutions, but is otherwise analytic. It permits an analytic construction 
of the microcanonical phase diagram of our system at small mass and charge. In the rest 
of this introduction, we will describe in detail our construction of small hairy black holes, 
their properties, and the phase diagram of our system. 

To begin with, we start our discussion with a consideration that may, at first, seem 
unrelated to the study of AdS black holes. Consider a spherically symmetric shell of a scalar 
field of frequency u> incident on a charged black hole in flat spacetime. One might naively 
expect a part of this wave to be absorbed by the black hole while the rest is reflected back 
to infinity. It is, however, a well known fact that the reflection coefficient for this process 
actually exceeds unity when u < e/j, (/i is the chemical potential of the black hole). Under 
these conditions more of the incident wave comes out than was sent in. This phenomenon, 
called superradiance [ll]], has immediate and well known implications for the stability of 
small RNAdS black holes, as we now explain. 

Consider a superradiant wave incident on a small charged black hole sitting at the 
centre of global AdS spacetime. Such a wave reflects off the black hole, propagates out to 
large r, but unlike the flat spacetime case, bounces back from the boundary of AdSs an d 
then finds itself re-incident on the black hole. This process continues indefinitely. As every 
reflection increases the amplitude of this wave by a fixed factor, this process constitutes an 



2 See Appendix H for a summary of notation employed in this paper. 
3 See however |7, I, B] for analytical studies in a related context. 

4 See lid] for earlier work on scalar condensation in black hole backgrounds in global AdSs. 
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instability of the charged black hole. A closely related instability, the so called black hole 
bomb, was discussed (in the context of a flat spacetime black hole surrounded by mirrors) 
as early as the 1970s @. 

As the spectrum of frequencies of a minimally coupled charged scalar field (in a gauge 
where A[ r ~°°^ = 0) in AdSs is bounded from below uj > Aq = 4 , we expect small charged 
black holes in AdSs space to exhibit superradiant instabilities whenever the condition 
e/i > to > Ao is satisfied 5 . Now the chemical potential fi of a small black hole is bounded 
from above by the chemical potential of the extremal black hole; i.e. fj? < fj? c = |. It 

follows that small charged AdS black holes are always stable when e 2 < ^ = e 2 = -y. 
When e 2 > e 2 , however, small black holes that are near enough to extremality suffer from 
a superradiant instability. 

The superradiant instability described above admits a very simple thermodynamical 
interpretation. Notice that the Boltzmann factor for a mode of energy Ao and charge e 
is given by e~ T 1 ( A °~ e ^ where T is the temperature of the black hole. Now this factor 
leads to an exponential enhancement (rather than the more usual suppression) whenever 
fie > Ao- In other words, a small charged black hole with fie > Ao is unstable against 
Bose condensation of the lightest scalar mode. Indeed the leading unstable mode of a small 
charged black hole with fi > — a is a small deformation of the lightest scalar mode in global 
AdSs space. 

The considerations outlined above suggests that superradiant tachyon condensation 
proceeds in the following manner. The black hole emits into a scalar condensate, thereby 
losing mass and charge itself. As the charge to mass ratio of the condensate (i.e superradiant 
mode), exceeds ^, the chemical potential of the black hole also decreases as this 
emission proceeds. Now the decay rate of the black hole is proportional to (Ao — fie) 
and so slows down as fi approaches It seems intuitively plausible that the system 
asymptotes to a configuration consisting of a fi « ^ stationary charged black hole core 
surrounded by a diffuse AdS scale charge condensate, i.e. a hairy black hole. We will 
provide substantial quantitative evidence for the correctness of this picture in this paper. 

In the discussion of the previous paragraph we have ignored both the backreaction of 
the scalar field on the geometry as well as the effect of the charged black hole core on the 
scalar condensate. However these effects turn out to be small whenever the starting black 
hole is small enough. In other words the end point of the superradiant instability of a small 
charged black hole is given approximately by a non-interacting mix of the black hole core 
and the condensate cloud at leading order. We will now pause to explain why this is the 
case. 

First note that the charge and energy density of the superradiant mode is contained in 
an AdS radius scale cloud. As the charge and mass the initial unstable black hole is small, 
the same is true charge and mass of the eventual the scalar condensate. Consequently, 
the scalar condensate is of low density and so backreacts only weakly on the geometry 

5 In Appendix |X| we verify this expectation by direct computation of the lowest quasi normal mode 
of this system. We find that for small R, the imaginary part of the frequency of this mode is given by 
3i? 3 (e/i — 4) = 3R 3 (e/i — Ao), where R is the Schwarzschild radius of the black hole. This imaginary part 
changes sign precisely where we expect the instability. 
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everywhere. For this reason the metric of the final solution is a small deformation of the 
RNAdS black hole with // = and the scalar condensate does not significantly affect the 
properties of the RNAdS black hole. On the other hand the condensate cloud is very large 
compared to the RNAdS black hole at its core. This difference in scales ensures that the 
charged black hole also does not significantly affect the properties of the scalar condensate. 

Motivated by these considerations, in this paper we construct the hairy black hole that 
marks the end point of the superradiant tachyon condensation process in a perturbative 
expansion around a small RNAdS black hole with \i = ^ and small but arbitrary radius. 
The perturbative procedure we employ in our construction is completely standard except 
for one twist, which we now explain. As is usual in perturbation theory, we expand out 
the metric, gauge field and scalar field in a power series in e which is the small parameter 
of our expansion 7 

i 2 

5V — 90fiu + e 92[mu + ■ ■ ■ 

A t = A ot + e 2 A 2t + ... (1.2) 

<j) = 601 + € 3 (j) 3 + . . . 

Here go^ u and Aq^ are the metric and gauge field of our starting RNAdS black hole so- 
lution. We then plug this expansion into the equations of motion, expand the latter in a 
power series in e, and attempt to solve the resultant equations recursively. Unfortunately, 
the linear ordinary differential equations that appear in this process do not appear to be 
analytically solvable in full generality. However it turns out to be easy to solve these equa- 
tions separately in two regimes: at large r (in an expansion in — which we call as the 
far-field expansion and mark by a superscript 'out') and at small r (in an expansion in 
r which we call the near-field expansion and mark by a superscript 'in'). Here r is the 
radial coordinate (that is zero at the black hole singularity and infinity at the boundary of 
AdS ) and R is the Schwarzschild radius of the unperturbed RNAdS black hole solution. 
The first expansion is valid when r 3> R, while the second expansion works when r<l. 
As we are interested in R <C 1, the validity domains of these two approximations overlap. 
Consequently, we are able to solve the resultant linear equations everywhere, in a power 
series expansion in R 2 . 

When the dust has settled we are thus able to solve for the hairy black holes only in 
a double expansion in e 2 and R 2 . This expansion is sufficient to understand small hairy 
black holes 8 . In section |2| below we have explicitly implemented this expansion to 0(e m R 2 ) 
for m < 5. Our calculations allow us to determine the microcanonical phase diagram of 

6 Note that, in contrast, for the small charged black hole at the core has its mass and charge concentrated 
within a small Schwarzschild radius. Consequently even a black hole of very small mass and charge is a 
large perturbation about the AdS vacuum at length scales comparable to its Schwarzschild radius. 

7 As we explain in Section ^ below, we find it convenient to choose e to be the coefficient of the ^ decay 
of the scalar field at infinity, i.e. the vacuum expectation value of the operator dual to the scalar field. 

s The technical obstruction to solving the equations at arbitrary R has a physical interpretation. Even 
at arbitrarily small e, the black hole reacts significantly on the condensate at finite R. Our system can be 
regarded as a non interacting mix of the black hole and the condensate only at very small R. Effectively, 
in this paper, we perturb around this non interacting limit. 
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our system, as a function of mass and charge at small values of these parameters ; our 
results are plotted for e = 5 in Fig.|] below (the results are qualitatively similar for every 
e provided e 2 > ttj an d ma Y a ^ so be simply generalised to the study of ( |hl| ) with a mass 
term added for the scalar field- see |6,4| ) . 




0.002 0.004 0.006 0.008 0.010 



Figure 1: Micro Canonical Phase Diagram at small mass and charge. The overlapping region is 
dominated by the hairy blackholc. 



As summarised in Fig.|T], hairy black holes exist only in the mass range 

(9e 2 - 64) „ , , N 3e / 32 \ 3 (5e 4 + 64e 2 - 1024) „ , , N 
Q+ K = 9 j Q 2 + 0(Q 3 ) < M < - 1 + — Q i — Lq2 + (Q 3 ), 



e lire 2 " " 16 \ 3e 2 y QAire 2 

(1.3) 

where M and Q are the mass and the charge of the hairy black hole. 

Above the upper bound in (|1.3| ) (i.e. in the shaded grey region above the blue line in 
Fig.|l|.), RNAdS black holes are stable and are the only known stationary solutions. The 
upper end of ( |1.3| ) coincides with the onset of superradiant instabilities for RNAdS black 
holes. The lower bound in (T3) is marked by the lowest (i.e. red) line in Fig.[l]. The 



9 The microcanonical ensemble is well suited to our purposes. We discuss the phase diagram in other 
ensembles, in particular the canonical and grand canonical ensemble, in Appendix [f] and ^ below. We are 
able to make less definite statements in these ensembles because it turns out that the system at a given 
fixed chemical potential and temperature often receives contributions both from small as well as big black 
holes. As the approximation techniques of this paper do not apply to big black holes, we are unable to 
quantitatively assess the relative importance of these saddle points. 
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extremality line for RNAdS black holes (the yellow line in FigJTj.) lies in the middle of 
this range in (1.3). At masses below lower bound of ( |1.3| ) (red line in Fig.H.), the system 



presumably has no states. 

As we have emphasised, a hairy black hole may approximately be thought of as a non 
interacting superposition of a RNAdS black hole and a scalar Bose condensate. At the 
upper bound of ( |1.3[ ) the condensate vanishes, and hairy black holes reduce to a RNAdS 
black holes. As we decrease the mass of a hairy black hole at fixed charge (or increase the 
charge at fixed mass), the fraction of the condensate increases. Eventually, at the lower 
end of ( |l.3| ) - the red line in Fig.|]. - the black hole shrinks to zero. In this limit 10 the 



solution reduces to a regular horizon free soliton (see [13[ for the appearance of a similar 
solitonic solution in a qualitatively similar context). The solitonic solution, on the red line 
of Fig.|l| is simpler than the hairy black hole solution. As R = on this solution, it may 
be generated as an expansion in a single parameter e; in Appendix ^ we have carried out 
this expansion to 17th order. Hairy black holes may be thought of as solitons with small 
RNAdS solutions in their centre (see the related general analysis 11 of [13] ). 

In summary, the hairy black hole interpolates between pure black hole and pure con- 
densate as we scan from the upper to the lower bound of Ql-3| ) (or down from the blue line 
to the red line in Fig. |].). Throughout the range of its existence, the hairy black hole is the 
only known stable solution 12 . It is also the thermodynamically dominant solution, as the 
entropy of the hairy black hole exceeds that of the RNAdS black hole of the same mass 
and charge, whenever both solutions exist. 

Note that the solitonic solutions described above have some similarities to so called 
boson stars 13 , which have been extensively discussed in the General Relativity literature, 
mainly in asymptotically flat space (see e.g. |l5| ) but also in asymptotically AdS spaces 
(see e.g. [jlq| ). However boson stars usually have scalar fields with a harmonic time 
dependence, which obstructs placing black holes in their centre (the scalar field would 
oscillate an infinite number of times as it approaches the horizon) . Our solitonic solutions 
are genuinely stationary 14 in a particular gauge. In this gauge the temporal component of 
the gauge field takes a particular non zero value at the origin of the soliton. This allows 
us to construct stationary hairy black hole solutions by placing charged black holes at the 



centre of the soliton (roughly via the procedure of [14].) if and only if the black holes are 



10 We take the R — > limit purely within classical relativity. Of course stringy and quantum gravity effects 
(including the one loop energy density of the gas outside the black hole) become important when the black 
hole becomes parametrically small. Such effects are important only in an infinitesimal wedge above the red 
line in Fig. [I] and depend on the detailed microphysics of the system (they would be different, for instance, 
in string theory and M theory). We ignore all such effects in this paper. We thank K. Papadodimas for a 
discussion on this issue. 

11 We thank G. Horowitz and H. Reall for drawing our attention to this reference. 

12 Except for excited solitons, see below for details. 

13 We thank G. Horowitz, M. Rangamani and K. Papadodimas for drawing our attention to the literature 
on boson stars and explaining their properties to us. K. Papadodimas has further drawn our attention to 
the fact that our soliton reduces precisely to a traditional boson star in the limit e — > 0. Of course hairy 
black holes exist only when e > e c and so do not exist in the small e limit. 

14 The presence of the gauge field allows us to evade Derrick's theorem. We thank K. Papadodimas 
pointing this out. 
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chosen so that their chemical potential matches the gauge field at the centre of the soliton. 
The last condition has a simple and intuitive thermodynamical interpretation; solitons and 
black holes can be put into an approximately non interacting mix only when their chemical 
potentials match! 

The hairy black hole we have described so far in this introduction is a weakly interacting 
mix of the RNAdS black hole and a condensate of the ground state of the scalar field. The 
reader may wonder whether it is possible to construct an excited hairy black hole solution 
that is a weakly interacting mix of a RNAdS black hole and an excited state of the scalar 
field. This is indeed the case. The set of spherically symmetric linearised excitations of 
a massless scalar field appear in a one parameter family 15 labelled by an integer n. The 
energy of the n th state is A n = 4 + 2n where n = . . . oo. It turns out to be possible to 
mimic the construction described above to construct excited hairy black holes that reduce, 
at small masses, to the superposition of a RNAdS black hole with [i = ^p- = 4 ~^ 2n with a 
condensate of the n th scalar excited state. It turns out that these excited hairy black holes 
are all unstable to the superradiant decay of the scalar mode with energy Ao = 4. They 
presumably decay to the ground state hairy black hole, in comparison to which they are 
all turn out to be entropically sub dominant. 

Each excited hairy black hole exists in a limited mass range which turns out to be a 
subset of the mass range (O). At the lower end of this range, each excited hairy black 
hole reduces to a horizon free scalar condensate, i.e. an excited soliton. Here we find a 
surprise. Recall that the instability of excited state hairy black holes is superradiant in 
nature. Superradiant instability rates scale like R 3 , and so go to zero as R — > 0. It thus 
appears that excited state scalars are actually stable to small fluctuations, and so cannot 
decay classically. The discussion of this paragraph actually leaves open the possibility 
that excited solitons may have an independent non super radiant instability mode. It is, 
however, possible to demonstrate that the spectrum of small oscillations about arbitrary 
unstable solitons has no exponentially growing eigenmode, at least within the assumption 
of spherical symmetry (see S.3). While this result does not rigorously prove the stability 
of excited solitons [17], |18| , it at least suggests that they are stable. We find this surprising 
and do not have a clear sense for its implications. 

It may be worth pausing to consider the relative merits and disadvantages of the per- 
turbative procedure of this paper versus the numerical approach more usually used to study 
hairy black branes. On the negative side our perturbative procedure gives us no information 
about the regime of large masses and charges (where the perturbative expansion breaks 
down). Within its regime of validity, however, our perturbative procedure is very powerful. 
It allows us, once and for all, to compute the phase diagram and thermodynamics of all 
relevant solutions - including each of the infinite number of excited state hairy black holes 
- as analytic function of the parameters of the problem (e.g. the mass and charge of the 
scalar field). Perhaps more importantly our procedure gives us qualitative intuitive insight 
into the nature of hairy black holes. For instance, as we have explained many times, the 
hairy black hole is an approximately non interacting mix of a RN-AdS black hole and the 

15 These modes are dual, under the state operator map, to the operator d 2n O, where O is the dimension 
4 operator that corresponds to the bulk field <f) according to the rules of the AdS / CFT dictionary. 
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scalar condensate. This picture together with a few lines of algebra, immediately yields a 
formula for the entropy of the hairy black hole, to leading order in its mass and charge (see 
|6.4| ). In other words the perturbative approach employed in this paper gives more than 
numerical answers; it helps us to understand why small hairy black holes behave the way 
they do. 

In this paper we have focussed on the black holes in global AdS at small mass and 
charge. Almost all previous studies of the Lagrangian Ql.l| ) have studied the system in the 
Poincare patch 16 . The local dynamics of black holes in global AdS reduces to the dynamics 
of black branes of Poincare AdS at large mass and charge. It follows that the phase diagram 
displayed in Fig.[l] should make contact with the results of previous analyses at large mass 
and charge. In section |7| we present a conjectural phase diagram that interpolates between 
the small mass and charge behaviour derived in this paper and the large mass and charge 
behaviour determined in previous work. The analysis of that section makes it clear that 
something special happens to the phase diagram of hairy at charges, in units of the inverse 
Newton constant, of order unity. We leave the detailed study of this special behaviour to 
future work. 

The reader who is interested in asymptotically AdS gravitational dynamics principally 
because of the AdS / CFT correspondence might legitimately complain that our choice of the 
Lagrangian fll.lj ) was arbitrary; the dynamics of charged scalar fields in any given example 



of the AdS/CFT correspondence is unlikely to be given by (IT). Our attitude to this is 
the following: we reg ard flOD as a toy model which we have chosen to study (in common 
with much earlier work on the subject), largely because it is a simple system that possesses 
several of the ingredients that are qualitatively important for hairy black hole dynamics. 
Our exhaustive study in this paper of the toy model sets the stage for a similar analysis of 
small hairy black holes in 'realistic' theories IIB theory on AdSs x 5 5 . We expect small of 
black holes in this special bulk theory to share many of the qualitative features discussed 
in this paper, but also to have some properties that result from dynamical features special 
to it. We have already commenced this study and hope to report on it in the near future. 

17 

Note added in v2 : After the first version of this paper appeared, we became aware 
of a recent work by Maeda.et.al.|2^] studying a similar subject. 



2. The basic setup for the Hairy BHs 
2.1 Basic equations of motion 

As mentioned in the introduction, in this paper we study the Lagrangian ( PI) . This action 
describes the interaction of a massless minimally coupled scalar field, of charge e, interacting 
with a negative cosmological constant Einstein Maxwell system. Through most of this 

16 See however for a study in global AdS. 

17 In this context it is however important to keep in mind that small black holes in IIB theory on AdSs x S s 
sometimes suffer from Gregory Laflamme instabilities (in addition to potential superradiant instabili- 
ties), an additional feature that complicates (but enriches) the dynamics of small black holes in a 'realistic' 
theory like IIB SUGRA on AdS 5 x S 5 . 
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paper we will be interested in stationary, spherically symmetric solutions of this system . 
However, in Appendix.^ we will generalise to the study of time-dependent configurations 
to investigate the stability against small fluctuations. We adopt a Schwarzschild like gauge 
and set 

ds 2 = -f(r)dt 2 + g(r)dr 2 + r 2 dQ 2 3 
A t = Air) 

(2-1) 

A r = Ai = 



The four unknown functions f(r), g(r), A{r) and (j)(r) are constrained by Einstein's equa- 
tions, the Maxwell equations and the minimally coupled scalar equations. It is possible to 
demonstrate that /, g, A, (j) are solutions to the equations of motion if and only if 

r (3/'(r) - 2e 2 rg(r)A(r) 2 (j)(r) 2 + rA'{rf) - 2/(r) ((6r 2 + 3) g(r) + r 2 (f)'(r) 2 - 3) = 
f(r) (3rg'(r) - 2g(r) (r 2 l <j>> \r) 2 + 3) 

+6 (2r 2 + 1) g(r) 2 ) - r 2 gir) (2e 2 g(r)A(r) 2 (f>(r) 2 + A'(r) 2 ) = 
rg(r)f'(r)A'(r) + f(r) (rg'(r)A'(r) + Ae 2 rg(r) 2 A(r)<j)(r) 2 

-2g(r) (rA"(r) + 3A'(r))) = 
g(r) ((r/'(r) + 6/(r)) <j>'{r) + 2rf{r)(j)"{r)) - rf(f)g'{r)<t>'{r) + 2e 2 rg(r) 2 A(r) 2 <j)(r) = 0. 

(2.2) 



The four equations listed in (2.2) are the rr and tt components of Einstein's equations, the 
Maxwell equation and the minimally coupled scalar equation, in that order. 

The equations (|2.2[) contain only first derivatives of / and g, but depend on derivatives 



upto the second order for <p and A. It follows that (|2.2j) admit a 6 parameter set of 



solutions. One of these solutions is empty AdSs space, given by f(r) = r~ + 1, g(r) 



2,i n („\ _ 1 



l+r 



7: 



A(r) = (f)(r) = 0. We are interested in those solutions to ( |2.2| ) that asymptote to AdS 
spacetime, i.e. solutions whose large r behaviour is given by 



f(r) = r 2 + l + 0(l/r 2 ) 
1 

l + r 2 ' ~ ' (2.3) 
Air) = 0(l) + 0(l/r 2 ) 

0(r) = 0(l/r A ) 



2(r) = ^2+0(l/r 6 ) 



It turns out that these conditions effectively impose two conditions on the solutions of 
Q2.2] ), so that the system of equations admits a four parameter set of asymptotically AdS 
solutions 18 . We will also be interested in solutions that are regular (in a suitable sense) in 
the interior. As we will see below, this requirement will cut down solution space to distinct 
classes of two parameter space of solutions (the parameters may be thought of as the mass 
and charge of the solutions). In particular, we are seeking the hairy black hole solutions 



3 For example, the equations above are easily solved in linearisation about AdSs ;the six dimensional 
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of the above equations that constitute the endpoint of the superradiant instability of small 
RNAdS black holes. To set the stage and notations for our computation we first briefly 
review the charged RNAdS black hole solutions in global AdS spacetime. 

2.2 RNAdS Black Holes and their superradiance 

The AdS-Reissner-Nordstrom black holes constitute a very well known two parameter set 
of solutions to the equations (|2,2|). These solutions are given by 



ds 2 
V(r) 

A(r) 
4>{r) 



-V(r)dt 2 + 



1 + r 2 - — 



dr 2 



+ r 2 dn 2 3 



1 + R 2 + -/x 2 



1 + r 2 + R 2 



2^R 



2 o R 4 

+ ~ 

3 

2l?2 



3 r 2 



(2.5) 



1? 



where [i is the chemical potential of the RNAdS black hole. The function V(r) in (£ 
vanishes at r = i? and consequently this solution has a horizon at r = R. In fact, it can 
be shown that R is the outer event horizon provided 



f<°-{l + 2R 2 ) 



(2.6) 



We will review later the thermodynamics of these solutions in more detail with a 
particular focus on small charged black holes whose R <C 1. Consider the small RNAdS 
black hole solutions of the system described by the Lagrangian in Ql.l| ). As we have 
explained in the introduction, in the limit i? < 1 we expect the solution in ( |2.5| ) to be 
unstable to superradiant decay provided e/x > Ao = 4. In appendix [A], we verify this 
intuitive expectation by determining the lowest quasinormal mode of this system in a 



power series in R. In a gauge where A t 
lowest mode is given by e~ lut where 



(r=R) 



0, we find that the time dependence of this 



u = (A - en) + R 2 (-6 + 3e/x - 4// 2 ) - i 3i? 3 (A - e/x) + 0(R 4 
= (4 - en) + R 2 (-6 + 3e^ - \\x 2 ) - i 3i? 3 (4 - e/x) + C(i? 4 ). 



(2.7) 



solution space is given by 



S/(r) = ai (l + r 2 )-^f 



«2 



r 2 (l + r 2 ) 2 

04 



50(r) 

The asymptotically AdS condition set ai = as — 



(2.4) 



r 

a 5 + as 



dr 



r 3 (l + r 2 
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Note in particular that 



Im (w) = -3i? 3 (A - efi) + 0(R 4 ) 



it follows that the time dependence e~ wt of this mode represents an exponential damping 
when fie < Aq but an exponential growth when fie > Aq. Consequently, small charged 
black holes are unstable when fie > Aq, in agreement with the intuitive expectations 
outlined in the introduction. Further, note that the decay (or growth) constant of the 
lowest quasi normal mode is given by 3i? 3 |A — e/x|, and goes to zero either when R goes 
to zero or as fi goes near As we have argued in the introduction, this motivates us to 
seek a hairy black hole solution which is constructed in a perturbation theory about these 
RNAdS blackholes. 

2.3 Setting up the perturbation theory 

The starting point of our construction is a small RNAdS black hole 



ds 2 = -V{r)dt 2 + 
V{r) = 
A(r) = fi 
at arbitrary but small R, and 



dr 2 



+ r 2 d^l 



2z?2l 



1 + r 2 + R 2 2 » R 



3 r 1 



(2i 



fi = fi{e,R) = J]e 2 > 2n (i?) 



n=0 
2k 



V2n(R) = ^2 H(2n,2k) R 
k=0 
4 

^(o,o) - - 



(2.9) 



Here fi = fi(R, e) is the as yet unknown chemical potential of our final solution. Note that, 
at the leading order in the perturbative expansion, fi = |. 
To proceed we simply expand every unknown function 



f(r,R,e) = J2^ n hn(r,R) 

oo 

g(r,R,e)=Y,£ 2n 92n(r,R) 

n=0 

oo 

A(r,R,e) = Y,t 2n A2n(r,R) 

n=0 

oo 

( f } (r,R,e)=Y J £ 2n+1 (f>2n + i(r,R) 



(2.10) 



n=0 
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Here fo, go and Aq are the values of the functions /, g and A for a RNAdS black hole with 
radius R and chemical potential fx = fio(R). given in $LQ). We expand our equations in 
a power series in e. At each order in e we have a set of linear differential equations (see 
below for the explicit form of the equations), which we solve subject to the requirements of 
the normalisability of 4>{r) and f(r) at infinity together with the regularity of 4>(r) and the 
metric at the horizon. These four physical requirements turn out to automatically imply 
that A(r = R) = i.e. the gauge field vanishes at the horizon, as we would expect of a 
stationary solution. These four physical requirements determine 4 of the six integration 
constants in the differential equation, yielding a two parameter set of solutions. We fix 
the remaining two integration constants by adopting the following conventions to label our 
solutions: we require that (j)(r) fall off at infinity like -% (definition of e), that the horizon 
area of our solution is 2ir 2 R? (definition of R). This procedure completely determines our 
solution as a function of R and e. We can then read of the value of \i in ( |2.9D on our 
solution from the value of the gauge field at infinity. 

As we have explained in the introduction, the linear differential equations that arise 
in perturbation theory are difficult to solve exactly, but are easily solved in a power series 
expansion in R, by matching near field and far field solutions. At every order in e we thus 
have a solution as an expansion in R. Our final solutions are, then presented in a double 
power series expansion in e and R. 

In the next few sections, we present a detailed description of the implementation of 
this perturbation expansion at order e and e 2 . In Appendix |B| we present explicit results 
for this perturbation expansion at higher orders. 



3. Perturbation theory at 0(e) 

We will now present a detailed description of the implementation of our perturbative expan- 
sion at 0(e). The procedure described in this subsection applies, with minor modifications, 
to the perturbative construction at 0(e 2m+1 ) for all m. 

In this section we wish to construct the first order correction around the black hole 

fo(r,R) = V(r), g (r,R) ' 



V(r) 



R 2 

A Q (r,R) =/x (l- -j) (3.1) 



V(r) = l + r 2 + ^ 



Plugging in ( p. 10 ), we expand the equations of motion in a power series in e to 0(e). Of 
course all equations are automatically obeyed at 0(e°). The only nontrivial equation at 
0(e) is D 2 (f) = where = V M — ieA^ is the linearised gauge covariantised Laplace 
equation about the background (|3.1|) . We will now solve this equation subject to the 
constraints of normalisability at infinity, regularity at the horizon, and the requirement 
that 4>(r) ~ -% at large r. 
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3.1 Far Field Region 

Let us first focus on the region r 3> R. In this region the background ( |3.1| ) is a small 
perturbation about global AdS space. For this reason we expand 

oo 
k=0 

where the superscript out emphasises that this expansion is good at large r. In the limit 
R — > 0, ( |3.lQ reduces to global AdS spacetime with At = -. A stationary linearised 
fluctuation about this background is gauge equivalent to a linearised fluctuation with time 
dependence e~ 4lt about global AdS space with At = (At is the temporal component of 
the gauge field). The required solution is simply the ground state excitation of a massless 
minimally coupled scalar field about global AdS 



I out 



v (i,o) - (i + r 2)2 
The overall normalisation of the mode is set by the requirement 

1 

r 

12 



r) = TT^m- (3-3) 



®)W = J J + 0(l/r B ). 



We now plug ( p,2j ) into the equations of motion D § = and expand to 0(R ) to 

out 
1,2' 



solve for (j)^- Here D 2 is the gauge covariant Laplacian about the background fl2.8|) . Now 



where (D 2 ) out is the gauge covariant Laplacian about global AdS spacetime with back- 
ground gauge field A t = |. It follows that, at 0(i? 2 ), 



1 



(l +r 2)2 



/ r\2\out tout / r\2\out tout I r\2\out 

K^O) V>(1,2) — ~\ L, l) 0(1,0) _ "l^J 

This equation is easily integrated and we find 
AOUt _2 (-3e 2 + 6 (e 2 - 32) (r 2 + l) log(r) - 3 (e 2 - 32) (r 2 + l) log (r 2 + l) - 32) 



3e 2 (r 2 + l) 3 
6e 2 - 64 \ /e (61og(r)r 2 - 3 log (r 2 + l) r 2 - l) 



+ M0.2 -3 



e d / V 6(r 3 + r) z 



(3.4) 



We could iterate this process to generate f^n^fe) till any desired fc. As in (|3.4|) , it 
turns out that the expressions ^"^fc) are mcreas i n gly singular as r — > 0. In fact it may be 
shown that the most singular piece of (ffflfizk) sca ^ es n ^e upto logarithmic corrections. 
In other words the expansion of <p out in powers of R 2 is really an expansion in -p- (upto log 
corrections) and breaks down at r ~ R. 
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In summary we have found that, to 0{R 2 ) 



6f*i 



r) 



1 

> 2 + l) 2 
"2 ' 



-3e 2 + 3fe 2 



32) (r 2 + l) log(r 2 /(r 2 + 1)) - 32) 



+ 



6e 2 - 64 



M(0,2) 

+ e>( J R 4 /r 4 ) 

The small r expansion of this result is given by 

" 4 



3e 2 ( r 2 + X )3 

e (61og(r)r 2 — 3 log (r 2 + l) r 2 



(3.5) 



6 (r 3 + r\ 



6f*(r) 



[1 _ 2r 2 + 0(r 4 )] + i? 2 
6e 2 - 64 



- 32) log(r) - 2 



1 + §)""- 



^(0,2) 



^2 + elog(r) + | + C(r 2 )] + 0(i? 4 ) 



(3.6) 



Note that this result depends on the as yet unknown parameter ^(0,2)- This quantity will 
be determined below by matching with the near field solution. 

3.2 Near Field Region 

Let us now turn to inner region r <C 1. Over these length scales the small black hole is far 
from a small perturbation about AdSs space. Instead the simplification in this region arises 
from the fact that background gauge field, which is of order unity, is negligible compared 
to the mass scale set by the horizon radius jjj. In other words the gauge field is a small 
perturbation about the black hole background in this region. To display this fact it is 
convenient to work in a rescaled radial coordinate y = and a rescaled time coordinate 
r = Note that the near field region consists of spacetime points with y of order unity. 
In these coordinates the background black hole solution takes the form 

2 i dy 2 2 r n 2 

+ y dQ 3 



ds z = R z ( -V{y)dr z + 
V{r) = 
A T = Rhq{1 - ^ 



1 



1 



V(y) 
3y z 



(3.7) 



1 



y- 



The explicit factor of R in A T in ( |3.7j ) demonstrates the effective weakness of the gauge 
field. This justifies an expansion of the near field solution in a power series in R 



(3.8) 



k=0 



To determine the unknown functions in this expansion, we must solve the equation D 2 (fi m = 
0, where D 2 is the gauge covariant Laplacian about the background ( |3.7D . Our solutions 



are subject to the constraint of regularity at the horizon. Further, they must match with 
the far field expansion in equations ( |3.5|) and (|3.6| ) above. 
Note 

(D 2 ) m = jp{Dl) in + {D 2 2 ) in + ... 

where (D 2 ) in in is the leading part of (D 2 ) m in an R expansion. At leading order we find 
Dq4> 1 q \(y) = The two linearly independent solutions of this equation are easily obtained 
by integration. The only solution regular at the horizon is the constant 

Co)(y) = 1 

where we have determined the value of the constant by matching with equations ( p.5[ ) and 
( |3,6| ) (more below about the matching). 

At next order in R 2 we obtain an equation of the form 

D 2 ^ 2) =-D 2 4>f 10) {y). 

Even though (f> l M Q )(y) is a constant, the RHS of this equation is nonzero because of the 
gauge field in (|3.1|) . The equation is easily solved by integration; imposing regularity of 
the solution at the horizon we find 19 

C 2) (y) =« + ^2 - 6e V - 128 lo s ( 3e2 - 32 ) lQ g { sf y ^ 32 ) ~ 192 log ^ 3g2y2 " 32 ) 



3e 2 (v 2 — 1) 

+ 6 log (3eV - 32) e 2 + 128 log ( J_ ^ ' ) log (3e 2 y 2 - 32) 

32 - 3eV 



+ 64 log 2 (3e 2 y 2 - 32) + 128Li 2 



32 - 3e 2 



(3.9) 



where Li2(x) is the polylog function as defined in Mathematica 6 



OO I. 

Z K 



Li2(^)=^T2 
fc=l 

The single unknown parameter a in this solution will be determined by matching below. 

The perturbative procedure described above may be iterated to arbitrary order. It 
turns out that the fields ^(i 2 m) a ^ high m are increasingly singular at large y. In fact it 
may be shown that the dominant growth of ^™ 2 m) 1S g ener i ca lly y 2m ■ It follows that the 
near field perturbative expansion is an expansion in (yR) 2 = r 2 . 



The apparent logarithmic singularities at y = 1, in two of the terms of (3.9), actually cancel 
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In summary 

<P\ n (y) =1 + R 2 a 
R 2 



+ 



3e 2 



6e 2 y 2 - 128 log (3e 2 - 32) log ( -^- % -L^ - 192 log (3e 2 y 2 - 32) 



6 log (3e 2 y 2 - 32) e 2 + 128 log j- ^^-Ji j log (3e 2 y 2 - 32) 



+ 64 log 2 (3e 2 y 2 - 32) + 128Li 2 



-32-3eV 



32 - 3e 2 



+ 0(yR) 4 



(3.10) 



The large y expansion of this result is given by 



bf{y) =l + R 2 



-2y 2 + -(e 2 -32)log(y) +a 



64vr 2 



+ 6 1 



32 



log(3) 



+ 3i2 ( - 192 lo § 2 (3333^2 ) + 3841og(3) log (3e 2 - 32) 



(3.11) 



+ 121og(e) (3e 2 + 64 log (3e 2 - 32) - 96) ) + O ( ^ 



In ( |3.11| ) we have determined <t> l i(y) in terms of the as yet unknown parameter a which 
will be determined by matching below. 

3.3 Matching 

In order to match the near and far field results, we substitute y = in ( |3.11| ) (the large 
y expansion of (ft 1 ™) and view the resultant expression as an expansion about small r and 
small R. As we have explained above, the resultant expression is reliable to all order in R 
but only to order C(r 2 ); all terms of order 0(r 4 ) or higher receive contributions from as 
yet undetermined fourth order terms in the perturbation expansion of (ft™ (3.11). 

We then compare this expression with the small r expansion of cp out , (p.6|) . We can 



generate this expansion to any order in r that we desire (merely by Taylor expanding 
( |3.6j )); however the resultant expression is clearly valid only to 0{R 2 ) in R (terms of 
0(R 4 ) obviously receive contributions from as yet undetermined fourth order terms in the 
expansion of cp out ). Terms of the form r°R°, r 2 R° and r°R 2 (together with logarithmic 
corrections) are reliably computed by both expansions and so must agree. The unknown 
parameters a and M(o,2) are determined to ensure this (as we have more conditions than 
variables we also obtain valuable consistency checks). We find 

6e 2 - 64 
M(o,2) = -3 



a 



2 (-9e 2 - 192 log (3e 2 - 32) + 288) 



9e 2 



log(3) 



2 (3e 2 - 32 log 2 (32 - 3e 2 ) + 32) 



3e 2 



+ 



64tH 
9?" 



18 (e 2 - 32) log(fl) + 61og(e) (3e 2 + 64 log (3e 2 - 32) - 96) 



This completes our determination of our solution to 0(R 2 ). The procedure described in 
this subsubsection can be iterated to obtain the solution at higher orders. 
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4. Perturbation theory at 0(e 2 ) 

We now briefly outline the procedure used to evaluate the solution at 0{e 2 ). We proceed 
in close imitation to the previous subsection. The main difference is that at this (and all 
even orders) in the e expansion, perturbation theory serves to determine the corrections to 
the functions /, g and A rather than the function (f>. The procedure described here applies, 
with minor modifications, to the perturbative construction at 0(e 2m ) for all m. 

4.1 Far Field Region 

In the far field region r>fiwe expand 



/r*(r) = £i? 2m /g m) w 

m=0 

oo 

&*(r) = Y. R2m 9™%m) (r) 

m=0 

oo 



(4.1) 



m=0 



Plugging this expansion into the equations of motion and expanding to 0(R 2m ) we find 
equations of the form 



i(r 2 {l + r 2 ) 2 9 ^ m){ r) 



d ( 7, 



out i 
(2,2m) 1 



2(1 + 2r 2 ) 



dr 



1 + r 2 



9& m) (r) + if 2>2m) (r) 



(4.2) 



— rv 



WhereP (2, 2m) ( r )' P (2,2 m) (ri 



~(2,2m)( r ) 

, (r) are the source terms at order e 2 R 2m which are 



P A 

r (2,2m)^ 

obtained from terms quadratic in <ft° ut . The most general normalisable solution to these 



equations takes the form 



5 (2,2m) 



r 2 (l + r 2 ) 2 r 2 (l + r 2 ) 2 



dx P (2,2m) ( X ) 



/(2!L)W = -(i + ^ 



dx 



-9 (2,2m) ( X ) + P (2,2m) ( X ) 



(4.3) 



A °(2,2m) ( r ) 



(2,2m) 



/•oo ^ r >>oo 
fc(2,2m) + / -3 / ^ P (2,2m)M 



Note that this solution has three undetermined integration constants b( 2 ,2 m ) > ^(2,2m) an d 



fc 



(2,2m)- 



Let us first focus on O(R ). The constants &(2,o) an( i ^(2,0) are determined by the 
requirement that the expansion of 5°2 0) an< ^ ^(20) a * small r starts out regular (i.e. has 
no term that goes like 4j). This requirement follows from matching with the near field 



- 18 - 



solution. For example, a term in g° 2 ^ °t results would match onto a term in gig 1 that 
scales like $ R2 . However g^(y,R) has a regular power series expansion in R and so does 
not have such a term 20 . The constant k( 2 ,2m) = ^{2,2m)i 1S as ve t undetermined. 

4.2 Near Field Region 

We now turn to the solution in the inner region r«l. As in the previous section we find 
it convenient to solve the equations here in the rescaled y and r coordinates. We expand 

oo 

f?(y) = J2 R2m fthrn)(y) 

m=0 

oo 

9T(y) = Y. R2m ^,2 m) (y) (4.4) 

m=0 

oo 

A^(y) = Y J R 2m ^t2 m) (y) 

The equations are slightly simpler when rewritten in terms of a new function 

fin / \ 

K {2 , 2m) (y) = V (y)g™ 2m) (y) + 

where 

v (y) = - A 

In terms of this function the final set of equations take the following form. 



dK {2t2m) (y) _ {K) 

dy - a {2,2m)\y> 

Ty [ V dy ) = S ^m) (f) + "(CO) [ Ty ) (4.5) 

Where S^^ m ^(y), S^\ m ^(y) an d S( 2 \ m )(y) are the source terms which depend on the 
solutions at all previous orders. The most general solution to these equations is given by 

dx S (2,2m)( X ) 



A V,2m)(y) =^(2,2m)(l ~ ^ + / ^ / dw ( 5 (£L») ( W ) + M(0,0) 4>) ( 

/fe™) (?/) =^f^^,2m) (y) + [ dx (S{i\ m) (x) + 2xK^ 2m) (x)) 



(4.6) 



20 At higher orders in the expansion in R, similar reasoning will not set 6(2, 2m) and /ip' 2m ) to zero but 
will instead determine them by matching with 3(2,2m-2V 
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In the above solution two of the four integration constants are chosen such that the solution 
obeys the requirement that A(l) = vanishes at the horizon (regularity of the gauge field 
at the horizon in Euclidean space) and that /(l) = (this is the requirement that the 
horizon is located at y = 1, which follows from our definition of R). It may also be shown 
that the remaining two constants in the inner solution at 0{R 2k ) may be determined by 
matching to the outer solution at the same order (0(R 2k )). 

In particular, the inner solution at order R° is completely determined by matching 
with the 0(R ) outer solution that we have already determined above, in terms of a single 
unknown M(2,o)- This yields the complete solution at order O(R ) in terms of this one 
unknown number. 

4.3 Iteration 

This process may now be iterated. Our determination of the inner solution at O(R ) 
permits an unambiguous determination of the integration constants in the outer solution 
at 0(R 2 ). This then allows the complete determination of the inner solution at 0(R 2 ), 
which, in turn, permits the determination of the outer solution at 0(R 4: ) and so on. This 
procedure may be iterated indefinitely. 

In summary the procedure described in this subsection permits the complete deter- 
mination of the 0(e 2 ) correction to our solution (order by order in R 2 ), as a function of 
the shift in the chemical potential \i at C(e 2 ), i.e. in terms of the as yet unknown num- 
bers jU(2,2m)- These numbers are left undetermined by 0(e 2 ) analysis, but turn out to be 
fixed by the requirement that there exist regular solutions of the scalar equation at C(e 3 ). 
This is completely analogous to the fact that the 0(e) shift in the chemical potential was 
determined from the analysis of the scalar equation at 0(e). 

It is relatively straightforward (though increasingly tedious) to carry out our pertur- 
bative expansion to higher orders in perturbation theory. The equations at odd orders in 
the e expansion serve to determine scalar field corrections, while equations at even orders 
serve to determine corrections to the metric and gauge field. In Appendix |b| we list explicit 
results for the correction to the metric, gauge field and scalar field at low orders in per- 
turbation theory. We will analyse the thermodynamics of these hairy black holes in more 
detail in later sections. 

5. The Soliton 

The hairy black hole solitons of the previous section appear in a two parameter family 
labelled by R and e. In general our solutions may be thought of as a small RNAdS black 
hole surrounded by a cloud of scalar condensate. As we have described in the introduction, 
the limit e — >• switches off the condensate cloud. In this limit (the blue line of Fig. [I]) 
hairy black holes reduce to RNAdS black holes. On the other hand, in the limit R — > the 
black hole at the centre of the condensate cloud shrinks to zero size, apparently leaving 
behind a horizon free scalar cloud. This is indeed the case. Indeed the solitonic solutions 
so obtained are considerably simpler than hairy black holes, as they may be generated as 
a single expansion in e. The linear differential equations that we encounter at every order 
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in this process are exactly solvable without recourse to the elaborate matching procedure 
described in the previous section. 

In this section we will directly construct the hairy black hole at R = in a perturbation 
expansion in e. We refer to the solution of this section as the 'soliton'. In order to construct 
the soliton, we search for all stationary charged solutions that are everywhere completely 
singularity (and horizon) free. We use global AdS as a starting point for these solutions, 
which we construct in a perturbative expansion in the scalar amplitude. As in the previous 
section we will only study rotationally invariant solutions, i.e. solutions that preserve the 
full SO (4) symmetry group of AdSs . 

At linear order the complete set of regular, asymptotically AdSs , 50(4) symmetric 
fluctuations about global AdS is given by 21 

" ,^,n 2F1 \-n, -(n + 2), 2, -r 2 ] , with uj n = 4 + 2n - ea 
n (l + r 2 )™+ 2 ( 51 ) 

A t = a, and Sg^ = 5A { = 0. 

The equation (|5.1[) is simply the most general rotationally invariant normalisable and 



regular solution to the equation d 2 (j). The constant a in (5A) can be set to any desired value 
by a choice of gauge. While ( jO] ) is generically time dependent, it reduces to a stationary 
solution when only one of the modes is turned on and o is chosen accordingly. i.e., 

2k + 4 

a = and 0^ oc o nk: 

e 

for some non-negative integer k. This yields a discrete set of stationary, one parameter, 
solutions to the equations of motion, labelled by their amplitude. In the rest of this section 
we describe the construction of the nonlinear counterparts of these solutions in a power 
series expansion in the amplitude. We refer to these stationary solutions as nonlinear 
'solitons'. 

Unlike the two parameter set of black hole solutions described in the previous section, 
the solitonic solutions of this subsection appear in a one parameter family, labelled by their 
charge; the soliton mass is a determined function of its charge. The fact that their are fewer 



21 We remind the reader that the Gauss's hypergeometric function 2-F1 [a, b, c, z] is a solution to the 
equation 

r / ri \ f rl \ / rl 

2-F1 [a, b, c, z]=0 



dz^ ) I dz ) \ dz ) dz 



defined by the series 



t-h r 7 , {a) k {b)k z k 

2F1 [a, b, c, z] = 2^ 



(c) k k\ 



where the 'Pochhammer symbol' (a)k is defined by the raising factorial 

(a) fe = a(a + l)(a + 2) . . . (a + k - 1) 

Note also that for n an integer, the function 2F1 [— n, — n — 2, 2, z] is actually an n-th degree polynomial in 
the variable z. 
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solitonic than black hole solutions is related to the fact that the solitons we construct in 
this subsection have no horizons, and therefore carry no macroscopic entropy. 

The ground state soliton (i.e. the soliton at n = 0) has a specially simple interpretation. 
It may be thought of as the nonlinear version of the Bose condensate, that forms when a 
macroscopic number of scalar photons each occupies the scalar ground state 'wave function'. 
It also represents the R — >• limit of the hairy black hole solution of the previous section. 
We will construct this solution in this section, postponing discussion of excited solitons to 
the next section. 

5.1 Perturbation theory for soliton 

To initiate the perturbative construction of the ground state soliton we set 

/(r) = l + r 2 + ^e 2n / 2 „(r) 

n 
oo 



9 ^ = TT72 +f2 e2n 92n( 
n=l 

A 00 

A{r) = - + Y j e 2n A 2n {r) 



n=l 

(5.2) 



6 n=l 



2\2 +^2<f>2n+i{r)e 



2n+l 



1 + 

v n=l 



and plug these expansions into (2.2). We then expand out and solve these equations order 
by order in e. All equations are automatically solved upto 0(e). At order e 2n the last 
equation in ( |2.2D is trivial while the first three take the form 
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On the other hand, at order e 2n+1 the first three equations in ( |2.2D is trivial while the last 
equation reduces to 

Here the source terms -P2r^( r )> ^2n( r )^ ^2n\ r ) an< ^ -^2n+i( r ) are the source terms which 
are completely determined by the solution to lower orders in perturbation theory, and so 
should be thought of as known functions, in terms of which we wish to determine the 
unknowns f 2n , 92n, Mn and (fon+l- 

Note that (||1|) are identical to the equations that appear in the far field expansion of 
the hairy black hole solution of the previous section. This is intuitive; in the limit R — > 
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all of the hairy black hole spacetime lies within the far field region. The soliton is simpler 
to construct than the hairy black hole precisely because it has no separate near field region. 
The differential equations that arise, at any given order of perturbation theory, may simply 
be solved once and for all, with no need for an elaborate matching procedure. 

The equations ( |5.3| ) are all easily integrated. It also turns out that all the integration 
constants in these equations are uniquely determined by the requirements of regularity, 
normalisability and our definition of e, as we now explain. 

The integration constant in the first equation of (|5.3| ) is determined by the requirement 
that g(r) is regular at the origin. The integration constant in the second equation is fixed 
by requirement of normalisability for /2n( r )- The constant for the first of the two integrals 
needed to solve the third equation is fixed by the regularity of ^2n( r ) at the origin. The 
constant in the second integral (an additive shift in A2 n ) is left unfixed at 0(e 2n ) but is 
fixed at 0{e 2n+1 ) (see below). 

The equation ( ]5.4| ) is also easily solved by integration. The constant in the first in- 
tegration needed to solve this equation is determined by the requirement of regularity of 
02n+i( r ) a t the origin. Once we have fixed this constant, it turns out that the solution for 
02n+i is generically non normalisable for every value of the second integration constant. 
In fact normalisability is achieved only when the previously undetermined constant shift 
of A<m takes a specific value, a condition that determines this quantity. The constant in 



the last integral that determines <fi from (5.4) is determined by our definition of e which 
implies that 

4> 2n +i ~ 0(l/r 6 ) 

for n > 1. 

5.2 Soliton upto C(e 4 ) 

In summary, the perturbative procedure outlined in this subsection is very easily imple- 
mented to arbitrary order in perturbation theory. In fact, by automating the procedure 
described above, we have implemented this perturbative series to 17th order in a Mathe- 
matica programme. 

We present some of the results, to this order, in Appendix [C]. In the rest of this 
subsection we content ourselves with a presentation of our results to 0(e 4 ). 
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(5.8) 

We will postpone the discussion of the thermodynamics of these solution for later. 
5.3 Excited state solitons and hairy BHs 

As explained around fl5.ip , the stationary solitonic solution constructed in the previous 
section is simply one (albeit a special one, as it has the smallest mass to charge ratio) of 
an infinite class of stationary solitonic solutions, each of which may be constructed in a 
perturbative expansion in e, exactly as in the previous section. We label solitonic solutions 
by an integer n; the n th stationary soliton has chemical potential fi = 4 + 2n at small 
amplitude. 

We have explicitly constructed the excited solitons with n = 1 and n = 2 upto a high 
order in perturbation theory. In Appendix [p] below we present some of the details of our 
results at low orders. Further, we can further construct a large class of excited hairy black 
holes which reduce to these excited state solitons as their horizon size goes to zero. These 
black holes may be thought of as a mixture of the excited solitons and a small RNAdS 
black holes with \i ~ +2 " . In Appendix.^, we construct this excited state hairy black 
hole at n = 1. We have a simple program in Mathematica that may be used to generate 
the excited hairy black hole solution at any fixed value of n. It should prove possible to 
generalise this construction once and for all at arbitrary n, but we have not attempted this 
generalisation. We will present a detailed analysis of the thermodynamics (and stability) 
of the excited solitons ad black holes later in the paper. 
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6. Thermodynamics in the Micro Canonical Ensemble 

In this section we compare the entropies of the various solutions constructed in this paper 
as a function of their mass and charge. We find it convenient to present all formulae in 
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terms of the rescaled mass m and the rescaled charge q. The physical mass and charge of 
the system, M and Q, differ from m and q by the rescaling 



M = — m 

8 



(6.1) 



The grand canonical partition function for the system is defined by the formula 

Z GC = Trexp [-T -1 {M - fiQ}] . 
where T is the temperature and /i is the chemical potential. 
6.1 RNAdS Black Hole 

The basic thermodynamics for an RNAdS blackhole is summarised by the following for- 
mulae 22 
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were Q is the charge, M is the mass of the black hole, S is its entropy, T its temperature 
and [i its chemical potential. We use the symbol Ah to denote the area of the outer horizon. 
The condition for R to be the outer horizon radius is 



^=^<\{l + 2R 2 ). 



(6.3) 



We are mainly interested in small RNAdS black holes with i?Cl. and the thermody- 
namic expressions can be simplified in this limit. The mass of RNAdS black holes at fixed 
charge is bounded from below; at small m and q, we have 

2 / ^- \ 3 

\<l) +0(q 4 ) 



m > 2\/-q + 



3 q , 



22 Throughout this paper, we find it convenient to consistently omit a factor of G 5 1 from all our extensive 
quantities. 
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For every pair (m, q) that obeys this inequality, there exists a unique black hole solution. 

At small mass and charge (with mass and charge taken to be of the same order) the 
entropy and the radius of the black hole is given by 

tt 2 R 3 



S 



m + \ m 2 — \q 2 

r2 = y_ f_ + 0(m 2 ; m(?) 



(6.4) 



At leading order in mass and charge, the chemical potentials of these black holes are given 
as solutions to the equation 

m If, 2/r 2 " 
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while the temperature is given by 
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where R? is given in (6.4). 



6.2 Soliton - Ground State and Excited states 

Using the soliton solution in the previous section, the mass of the ground state soliton can 
be easily determined as a function of its charge. We find 



rn 



16q 2 / 64\ 2 . 



(6.7) 



In Appendix C.l, we give the relation of the mass and the charge implicitly upto higher 
orders. 

Upon continuing to Euclidean space, our soliton yields a regular solution for arbitrary 
periodicity of the Euclidean time coordinate (it is similar to global AdS spacetime in this 
respect). It follows that, within the classical gravity approximation, this soliton can be in 
thermodynamical equilibrium at arbitrary temperature. As the soliton has no horizon, its 
entropy vanishes in the classical gravity approximation. This implies that the free energy 
of the soliton is equal to its mass. 

The chemical potential of the soliton is given by the value of the gauge potential at 
infinity 

A /n RA \ 

+ 0(q 2 ) (6.8) 



+ 



64 
7? 



Note that the coefficient of q in the formula above is positive when e 2 > ^ = e 2 so that the 
chemical potential of the soliton increases with charge whenever hairy black holes exist. It 
is plausible that the only classical gravity state in the system with [i < - is the vacuum 
(or more precisely a thermal gas about the vacuum; this gas is absent in classical gravity). 
The grand free energy of the ground state soliton is given by 

3437re 2 (e/u-4) 4 



G(p) = M -TS - nQ 



4 (9e 2 - 64)' 



+ 0((^-4/e) 3 ) 



(6.9) 
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The above analysis is easily generalised to excited state solitons. For the general excited 
state solitons, we present thermodynamical formulae only at leading order. The mass and 
chemical potential of the soliton are given by 

m _ 4(4 + 2^ 

4 4- 1 „, < 6 ' 10 ) 
f = — — + 0(g) 
e 

We have explicitly constructed the excited solitons with n = 1 and n = 2 upto a high 
order in perturbation theory. In Appendix [D| we present some of the details of our results 
at low orders. This allows us to give the thermodynamic formulae for these cases upto a 
higher order At n = 1 we find 

3"9 , 1 2544 \ _2 , „, 3 , 
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Note that the coefficient of q in the expansion of \i in the expansion of ( |6.11| ) is positive 
whenever e 2 > 24 so that the first excited hairy black hole exists. On the other hand, the 
coefficient of q in the expansion of /x is negative at e 2 = , the threshold for the existence 
of the second excited hairy black hole (see below). 

6.3 Dynamical Stability of solitons 

In this subsection, we comment on the dynamical stability of the excited solitons. In par- 
ticular we will prove below that the spectrum of small fluctuations about excited solitons 
have no 50(4) symmetric exponentially growing modes, within the e perturbative expan- 
sion. This result suggests (but does not strictly prove |l7], [l8|]) that small excited state 
solitons are all dynamically stable against small fluctuations. 

As we have mentioned above, the normal modes of the scalar field constitute the only 



50(4) symmetric fluctuations of (LI) about global AdS spacetime. At small e the solitonic 
solution is everywhere a small perturbation around global AdS spacetime. It follows that 
the 50(4) symmetric perturbations about the soliton, at small e, are small perturbations 
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of spherically symmetric scalar normal modes about global AdS spacetime. These modes 
obey the equation 

D 2 (p = (6.13) 

where D is the gauge covariant derivative about the soliton background. We study pertur- 
bations of the form 

<f>(r,t) = ^(r)e -<w * 

and wish to investigate whether the frequencies u (which are all real about global AdS ) can 
develop a small imaginary piece about the solitonic background. We will now demonstrate 
that this is impossible in the e expansion. To establish this we multiply the equation ( 6.13|) 
by 4>* and integrate the resultant scalar over AdS spacetime. We find 

J ^g\g 00 \(u-eA t (r)) 2 \4,\ 2 = J ^gg rr \d r ^\ 2 

Here g^ u is the soliton metric and A t {r) is the gauge field of the solitonic solution. 
Now recall that A t = ± + 0(e 2 ). It follows that 



j^gg rr \d^\ 



2 



= v^_n , :, +0(0- 



As the leading term on the RHS is 0{e ) and positive, it follows that u is real within the e 
expansion. Consequently the spectrum of spherically symmetric small fluctuations about 
the soliton background does not have exponentially growing modes in the e expansion. 
This suggests that all excited solitons are classically stable. We find this result surprising, 
and think that it warrants further study. 

6.4 Massive scalar : Hairy black hole thermodynamics 

We concluded at the end of the previous subsection that the excited solitons seem to 
be classically stable. In contrast, our calculations in the appendix |A| indicate that the 
RNAdS blackhole can become super radiantly unstable. It is an interesting question to ask 
what is the thermodynamics of the resultant hairy black hole. Since, we have an explicit 
construction of the hairy black hole solutions, we can directly go ahead to compute the 
thermodynamic quantities for these hairy solutions. Before doing that however we wish 
to present an argument in this subsection which gives us some intuition about the kind of 
thermodynamics we should expect at the leading order. 

We will present this argument in a slightly more general framework than we have 
been working till now - we wish to consider the effect of adding a scalar mass term to the 
Lagr angian ( |1 . 1|) . , i . e . , we work with a more general system 



\ (K[g] + 12) - \T^ V - \D^\ 2 - m\\$ 



.14) 



Ffiv = V^A/ - V V ^ M and = V M - ieA^ 



- 28 - 



This system has a minimally coupled charge scalar with mass ms and charge e in AdSs. By 
the standard rules of AdS/CFT, the dual boundary operator has a scaling dimension 




2 + JA + ml 



d=4 



In a gauge where A r t = °° = 0, this is also the energy of the lowest <j> mode in vacuum AdSs. 
For the case = 0, this reduces to Ao = 4. 

The other spherically symmetric modes of <p (dual to the descendants d 2n O§ )have an 
energy 

A n = A + 2n = 2 + + + 2n 

For the case = 0, A n = 4 + 2n. Hence, in a gauge where A r t =co = 0, the energy of 
the n-th excited state is also given by A n . We can form a non- linear Bose condensate 
by dumping a charge Q so i into this n-th excited state - this is equivalent to populating 
this excited state with Q so i/ e number of charged Bosons. To the leading order, where we 
neglect self-interaction between these Bosons, the mass of such a soliton is given by 

M sol = ^A n + 0{Ql ol ) 



This sets the chemical potential of the soliton to be 
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The entropy of such a solution is zero S so i = 0. This in particular means that within this 
approximation, this solution exists at arbitrary temperatures T so [. 

Now, let us form a hairy black-hole by placing at the core of this non-linear Bose 
condensate a small ordinary RNAdS black hole with a small outer horizon radius R and 
chemical potential hbh- Such a black hole has a mass 
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an entropy 



and a temperature 
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If the number of Bosons Q ao ije is small, then the condensate outside is a small per- 
turbation on the RNAdS black hole. And if the radius R of the blackhole is small, then 
it is a small perturbation on the soliton on length scales large compared to R. Hence, if 
both these conditions are met, it is legitimate at the leading order to assume that there is 
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no interaction between the core and the condensate parts of the hairy black hole. In this 
regime, since the core and the condensate can still exchange charge and energy, all that is 
needed for a stationary solution is that the core and the condensate be at a thermal and 
chemical equilibrium,!. e., 



T = T so i = Tbh 
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and 
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Using these equilibrium conditions we want to figure out the 'mole fractions' of these 
two phases at equilibrium as a function of total mass and charge 



M = M so i + M BH and Q = Q sot + Q BH 



This is easily done and we get the mass fractions of the core and the condensate are given 
by 
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The charge fractions of the core and the condensate are given by 
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The radius of the black hole at the core is 
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and the entropy of the hairy black hole is given by 
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The existence region of the hairy black holes is in between where the hairy black hole 
coincides with the RNAdS black hole on one side and where it coincides with the pure 
soliton on the other side. This gives the existence region as 
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where e c is the critical charge above which the pure black hole becomes superradiantly 
unstable to radiation in the n-th excited state. 

In this regime, the upper bound on M is a decreasing function of A„ whereas the 
lower bound is an increasing function of A n . This implies that the existence region of n-th 
excited state hairy black hole is entirely inside the existence region of (n-l)th excited state 
black hole (see Figure]2|). Further, in this regime, one can show that the radius R is a 
decreasing function of A n . Hence, the higher excited state hairy black holes have smaller 
cores and consequently are entropically subdominant to the lower excited state hairy black 
holes. 

Armed with the above intuition, in the next few sections, we will derive the thermo- 
dynamics of hairy black holes with = directly from our solutions and show that their 
leading order behaviour is captured by the kind of non-interaction arguments that we have 
presented in this subsection. 

6.5 Ground State Hairy Black Hole 

Once we have our solutions for hairy black holes from Appendix [b], the evaluation of their 
thermodynamic charges and potentials is a straight forward exercise. At low orders in the 
perturbative expansion we find 23 
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23 Throughout this paper, we find it convenient to consistently omit a factor of G 5 1 from all our extensive 
quantities. 
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It may be verified that these quantities obey the first law of thermodynamics 



dM = TdS + fidQ. 



Equation ( |6.20 ) above lists formulae for the mass and charge of small hairy black holes as 
a function of R and e. Inverting these relations we find 
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Hairy black holes exist for all positive values of R and e. Of course R 2 and e 2 are 
positive; this implies that the mass and charge of hairy black holes vary over the range 
fll.3| ). As we have mentioned in the introduction, it is possible to satisfy this inequality 

only when e > = e c . Assuming this is the case, we have hairy black hole solutions 
only within the band ( |1.3| ), At the upper end of the band the solution reduces to a RNAdS 
black hole (e = 0). At the lower end of the band the solution reduces to the soliton with 
R = 0. 

It is now a simple matter to plug ( 6.21| ) into ( |6.20| ) to determine the entropy, temper- 
ature and chemical potential of the black hole as a function of its mass and charge. We 
find 
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As we have explained in the previous subsection, at leading order, these formulae have 
a very simple and intuitive explanation in terms of a noninteracting mixture of a RNAdS 
black hole and the ground state soliton. It is easily checked that the formulae in this section 
agree with the expressions that we had derived before if we put A n = 4. 
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6.6 Excited Hairy black holes 

We begin by reporting some of the basic thermodynamical formulae that follow from the 
formulae presented in Appendix.^. For the mass, charge and chemical potential we find 




Figure 2: Microcanonical ensemble for e = 5.4 : Existence region of ground state and excited 
state hairy black holes. 

The n th excited state is a small deformation of a small RNAdS black hole at /x = 
4 + 2 " . RNAdS black holes at this value of the chemical potential have n — 1 superradiant 
instabilities. It follows that the n th excited hairy black hole also has n — 1 unstable linear 
fluctuation modes, which tend to flow the black hole to lower excited (generically ground 
state) hair black holes. 

The n th excited hairy black hole exists only when e 2 > 2 ( 4 + w ) _ When this condition 
is fulfilled, the n th excited state black hole exists only when 

8, , /9 , (3e 2 + 8(n + 2) 2 ) g , 9n 

a<" + 2 >" + °<«"> s £ m„ + 2) + 0<q) 

This is completely in accordance with our non-interaction argument as expected. 
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7. Discussion 



In this paper we have demonstrated that very small charged hairy black holes of the 
Lagrangian (|0|) are extremely simple objects. To leading order in an expansion of the 
mass and charge, these objects may be thought of as an non interacting superposition of a 
small RNAdS black hole and a charged soliton. The different components of this mixture 
interact only weakly for two related reasons. The black hole does not affect the soliton 
because it is parametrically smaller than the soliton. The soliton does not backreact on 
the black hole because its energy density is parametrically small. 

We have constructed an infinite class of hairy black hole solutions labelled by a single 
parameter n. Excepting the ground state hairy black hole, each of these solutions is 
classically unstable. The time scale for this instability is proportional to the area of the 
RNAdS black hole that sits inside the hairy solution, and goes to zero in the limit that this 
area goes to zero. In particular all excited state solitons could well be stable configurations 



(see 3.3). We find the likely existence of an infinite number of classically stable classical 
solutions surprising, and do not have a good feeling for the implications of this observation. 
Of course these excited solitons will all eventually decay to the ground state hairy black 
hole via quantum tunnelling, but the rate for this decay will be exponentially suppressed. 
It would be interesting to construct the instanton that mediates this decay process. 

We have constructed hairy black holes in a perturbative expansion in their mass and 
charge. As we increase the mass and charge, the soliton and the black hole begin to 
interact with each other. At large mass and charge (where this system has been intensively 
previously investigated) there is probably no sense in which the hairy black hole can usefully 
be regarded as a mix of two independent entities. In fact we suspect that the soliton does 
not even exist as an independent object at large enough charge. It is very natural to 
wonder how the phase diagram of Figure 1 continues to large mass and charge. We sketch 
one possibility for this continuation in Fig. |3| 




Figure 3: Proposed Microcanonical phase diagram for large M and Q. 



The main point of interest of the conjectured phase diagram of Fig. ||| is the lower 
edge of the diagram. As we have explained below, the hairy black hole phase is bounded 
from below by the solitonic solution at small mass and charge. The temperature of the 
hairy black hole also tends to infinity as we approach this line. On the other hand, it 
seems plausible that the solitonic solution goes singular past its 'Chandrasekhar' critical 
charge q = q c (the dot in Fig. |3|). If this is indeed the case, it is of great interest to know 
the nature of the lower bound of Fig. || at higher values of the charge (the jagged line in 
Fig. H). We get some information from studies of the system in the Poincare patch limit 



(see for instance [21, |22j, |23j),so asymptotically large charges in Fig. y. In this limit prior 
studies appear to suggest that the jagged line is the zero temperature limit of hairy black 
brane solution 24 . If the phase diagram of Fig. || is indeed correct, it is natural to suppose 
that the jagged line everywhere represents the zero temperature limit of a hairy black hole. 
This suggests that the neighbourhood of the limiting soliton solution (the dot in Fig. ||]) 
is extremely interesting. Hairy black holes on the solid lower line in the neighbourhood of 
this dot are at infinite temperature. On the other hand, points on the jagged line in the 
neighbourhood of this dot are at zero temperature. As we circle the dot from the solid to 
the jagged line we presumably pass through all temperatures in between. While all this is 
very speculative, would be very interesting to investigate it further. 

Although all the analysis of this paper focussed on the concrete and especially simple 



context of charged black holes in the system (1.1), the basic physical picture of small hairy 



black holes as a linear combination of approximately non interacting pieces is based on 
very general considerations, should apply equally to the study of any bulk gravitational 
asymptotically AdS system that hosts RNAdS black holes which suffer from a superradiant 
instability. It should be straightforward to generalise the calculations of this paper to 
systems in which ( |l.l| ) is modified by the addition of a potential for the scalar field, and / 
or is studied in different dimensions. It may also be possible to generalise the constructions 



of this paper to superradiant instabilities in charged rotating black holes |24|, P5[ . 

As we have explained in section at leading order in perturbation theory, small hairy 
black holes exist for e 2 > but do not exist when e 2 < The case e 2 = ^ lies on 
the edge, and is particularly interesting. The question of whether hairy black holes exist 
at this critical value is determined by a second rather than first order calculation. In the 
system under study in this paper, it turns out that hairy black holes do exist for e 2 < -y. 
One way of understanding this statement goes as follows. 

At any fixed value of the charge q, hairy black holes exist if and only if e 2 > e 2 (g) 
where 

2/ \ 32 16 „, 2n 
el{q) = ^-^q + 0{q 2 ) 



24 It is interesting that the zero temperature limit of hairy black branes appears to depend qualitatively 
on the mass of the charged scalar field, and m — case is rather special. In the small black hole limit, on the 
other hand, we do not expect a qualitative dependence of our phase diagram on the mass of the scalar field. 



Indeed the leading order thermodynamics of hairy black holes in the Lagrangian (1.1) supplemented by a 
mass term for the scalar field can be easily obtained based on general non-interaction arguments (see 3.4) 
and it gives results that are qualitatively similar to the m — case. We thank A. Yarom for a discussion 
about this point. 
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(one may derive this result by comparing the mass of the soliton with the extremal RNAdS 
black hole at equal charge). It follows that small charge hairy black holes do exist at 
e 2 = tP, but we have to go beyond leading order in perturbation theory to see this 25 . 
Although we will not elaborate in detail in this paper, it turns out that the properties of 
small charged black holes at e 2 = ^ differ qualitatively from the properties of the same 
objects at larger e 2 . This observation is particularly relevant as it turns out that charged 
scalar fields are forced to sit at this critical value of charge in some natural supersymmetric 
bulk theories. We postpone the elaboration of these remarks and their consequences to 
future work. 

A very interesting arena in which the ideas of this paper might find application is 
to the study of charged and rotating black holes in IIB theory on AdSs x S 5 . Small 
charged rotating black holes in IIB theory on AdSs x 5" 5 sometimes suffer from superradiant 
instabilities. It should be possible to apply the ideas of this to the study small black holes 
on AdSs x S 5 (we have started work in this direction). In this context one should, however, 
also keep in mind that small black holes in this system sometimes also suffer from Gregory 
Laflamme type instabilities [19] even when they are uncharged, a feature that complicates 



(but enriches) the study of this system. 

Recall also that IIB supergravity on AdSs x S 5 hosts a 4 parameter set of one sixteenth 



BPS supersymmetric black hole solutions^, 27, 28]. It is possible that there exist new 
BPS hairy black holes consisting of a non interacting mix of these SUSY black holes with 
a SUSY graviton condensate. It would be very interesting to investigate this further. 

Acknowledgements 

We would like to thank R. Gopakumar, S. Hartnoll, H. Liu, G. Mandal, K. Papadodi- 
mas, W. Song, A. Strominger, T. Takayanagi, S. Trivedi, S. Wadia and A. Yarom for 
useful discussions. We would also like to thank S. Hartnoll, G.Horowitz, K. Papadodimas, 
M.Rangamani, H.Reall , A. Strominger, S. Trivedi and A. Yarom for useful comments on a 
draft version of this manuscript. SM would like to thank IPMU for hospitality while this 
work was being completed. The work of SM was supported in part by a Swarnajayanti 
Fellowship. We must also acknowledge our debt to the steady and generous support of the 
people of India for research in basic science. 



A. Superradiant Instability of small black holes 

In this section we analyse the dynamical stability of RNAdS black holes in AdSs to super- 

25 On the other hand, the study of the near horizon BF bound in highly charged extremal black hole 
backgrounds (which are locally well-approximated by extremal black-branes) indicates that as q — > oo, 

Note that the leading deviation from the black-brane result of e 2 = 3 is positive and in the large q limit e 2 
continues to be monotonically decreasing. In other words, all available data is consistent with the conjecture 
that e^.(q) is a monotonically decreasing function that interpolates between 4^ and 3 as q ranges from to 

CO. 
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radiant emission in the presence of a massless charged minimally coupled scalar field. As 
we have explained in the introduction, we intuitively expect very small black hole to be 
unstable whenever /i > ~. In this Appendix we verify this expectation by computing the 
frequency of the lowest quasi normal mode of the black hole in a perturbative expansion 
R, the radius of the black hole. We find that the imaginary part of this frequency flips 
sign (from stable to unstable) as \i increases past |, exactly as we expected on intuitive 
grounds. 

We wish to compute the lowest quasi normal mode of (|2.8| ) at small R. By definition, 
quasi normal modes are regular at the future horizon, so it is useful to work in coordinates 
that are good at the future horizon. We choose to work in ingoing Eddington Finkelstein 
coordinate; in other words we replace the Schwarzschild time t with the new time variable 



where V(r) given by 



r 2 J y 3 f 
In these new coordinates the background ( |2.8| ) takes the form 

ds 2 = 2dvdr - V(r)dv 2 + r 2 (d0 2 + sm 2 (8)dip 2 + sin 2 (6l) sin 2 (tp)dX 2 ) . (A.2) 

with 

R 2 \ 

(A.3) 



t+ 1 W) dr ' 



V (r)= (1-3 (l + r 2 + i? 2 -^ 2 ) (A.l) 



A v = p ( 1 r2 . . 



A —JL-U-E 

" V (r) V r 2 

A linearised scalar fluctuation about this background takes the form 

D^D^(v,r) = Q, (A.4) 

where 

£> M = V M - ieA^, 



with being the background gauge field ( A.3| ) 



In the rest of this appendix we will solve (A.4) separately in a far field region, r>i? 



and a near field region r<l. In the limit R <C 1, of interest here, the solution may then 
be determined everwhere by matching the two solutions in their overlapping domain of 
validity. The matching procedure may be carried out systematically in a power series in 
R, and turns out to determine the frequency of the quasi normal mode in a power series 
in R. We carry out this procedure to order R 3 , the first order at which the quasi normal 
frequency develops an imaginary component. 26 



26 The fact that the quasi normal frequency first develops an imaginary piece at 0(R 3 ) is simply related 
to the fact that the area - and so low frequency absorption cross section of a black hole in 5 dimensions - 
scales like R 3 . 
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A.l Solution in near field region 

When r <C 1 it is useful to work with the rescaled coordinate y given by 

r = Ry. 

Let the scalar fluctuation take the form 

$ in {v, y) = exp (-iojv) <l> in (y). (A.5) 

where 

<£ m (y) = % n (y) + <S>\ n (y)R + $ 2 n (y)R 2 + 0(R 3 ). (A.6) 

and 

w = 4 - + + w^i? 2 + u/ 3 ^ 3 + 0(i? 4 ) (A.7) 

Note we have chosen to study the quasi normal mode with the frequency 4— fie+0(R); here 
we have used the physical expectation that the lowest quasinormal mode should reduce to 
the lowest normal mode in the limit R — >• 0. The energy of the lowest normal mode is 4. 

It is a simple matter to solve ( |A4| ) perturbatively in R. Imposing the physical require- 
ment of regularity at the horizon we find 

% n (y) = d , 

mv) = 6(2 ^_ 3) ( 6 ( 2 ^ 2 - 3 ) ^ - - 4 )) 



+»do(e/z - 4) (log (3y 2 - 2/x 2 ) - 2 log(y + 1)) + 4\/6/u 3 tanh" 1 

&2 n (y) = -\d y 2 (eV-8e/x + 20) 

y(4 - efi) (do(e/i - 4) (2\/67r^ 3 - 91og(3)) + 6idi (2/x 2 - 3)) 
+ 6(2/i 2 -3) 

- 4d (en - 2n 2 - 3) log (-) + d 2 + o(- 



n 



(A.8) 

where do, &\ and d 2 are as yet undetermined integration constants (they will be determined 
below by matching). For brevity we have also presented the result assuming uft^ = 0, a 
result that turns out to be forced on us by matching with the far field expansion below. 

A. 2 Solution in the far field region 

In the outer region the fluctuation takes the form 

r ) = exp (-w) <S> out (r). (A.9) 

where 

$ oui (r) = <S>° ut (r) + $f*(r)il! + <S>° 2 u \r)R 2 + § z u \r)R z + 0(R 4 ). (A.10) 
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Solving the equation of motion subject to the requirement of normalisability at large r we 
find 

p— i(efi— 4) tan - 1 (r) 

$J"*(r) = 0, 

-i(e^-4)tan _1 (r) , o 

^ (r) = — - -3 3ir 2 (2/x 2 + 3) (e/i - 4) + - (r 2 + l) r (-8e/x log (r 2 + l) 
or (r z + 1) v <s 

+16 (e/i - 2/l 2 - 3) log(r) + 2i (e/i (2/i 2 + 9) - 8 (2/x 2 + 3)) tan _1 (r) - 2ivre/i 3 

-9z7re/i + 16/i 2 log (r 2 + l) + 24 log (r 2 + l) + lGivr^ 2 + 24ivr) 

+2i (2/l 2 + 3) (e/i - 4) - 4r (2/i 2 + 3) 



(A.11) 



and, 



g— i(efi— 4) tan 1 (r) , 

= g — w 3 (3i7rr 2 + 3r 2 (log (r 2 + l) - 21og(r) - 2itan~ 1 (r)) + 1 

6 (r 3 + r) ^ 

(A.12) 

where once again we have presented the results only for u)W = 0. We have also plugged 
a/ 2 ) = —6 + 3e/i — 4/i 2 in the expression for ^(V) and $3(r) (this is forced on us by the 
matching condition below). Here besides imposing normalisability at infinity, we have also 
demanded that the coefficient of the leading normalisable piece is one. 

A. 3 Conditions for patch up 

In order to complete our determination of the solution, we must now match the near and far 
field solutions. The logic for this matching procedure is exactly as described in subsection 
|3.3| . Implementing this procedure we find 

d = 1, 

i(e/i - 4) (2s/6iTTfj, 3 - 91og(3)) 



d 



6 (2/x 2 -3) 



d2 = 12 ( 4 ^ + 3 ) ~ ^ + M ^ + 48 ^ ~ 2 ^ ~ ^ l ° g ^ 



(A.13) 

+ 14) + 48 (e/i - 2/i 2 - 3) \og(R) 
-3iir (2e/i 3 + 9e/i - 16/i 2 - 24) V 

Also the quasi-normal frequency is determined to be 

u = 4 - e/i - R 2 (6 - 3e/i + 4/i 2 ) - R 3 (3t(4 - e/i)) + O (R A ) . (A.14) 

Once these matching conditions are imposed, the large y expansion of the near field solution 
(with y substituted by -^) and the small r expansion of the far field solution both share 
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the common expansions 



<5> ou \r) = (l - ir(en - 4) + r 2 (~ t 



e 2 n 2 + 4e/i - 10 ) + O [r 



( i (2/j, 2 + 3) (e/z - 4) , , , , . 

+ ^ " + (4(//(e - 2 M ) - 3) log(r) 



+ -(2// 2 + 3) (e/i(e/i - 8) + 14) 



i»7r (2e/i 3 + 9e^ - 16^ 2 -24)^+0 {r)\ R 2 

fezi) +0 fI)W +0(R4) . 



$ m (r) = (r 2 I - - e V + 4efi - 10 j - ir{e^ - 4) + 1 + O 

+ ( ( 4(>(e - 2a*) - 3) log(r) + i (2,u 2 + 3) (e/i(e/i - 8) + 14) 



7 ivr (2e/i 3 + 9e/u - 16^ - 24) ) + ' v r ~ ' v ~ r ~ - + O [ 



(2 M 2 + 3) (e/i-4) + / 1 

V 7" 



i{en - 4) 
27 1 



R 3 + (R 4 ) . 



(A.15) 



i? 2 



Equation (A. 14) is the main result of this Appendix. Note that the imaginary part of 
oj turns positive as \i exceeds |, demonstrating that RNAdS black holes with fi > | suffer 
from a super radiant instability. 



B. Results of the Low Order Perturbative Expansion of the Hairy Black 
Hole 

In this Appendix we present explicit results for the perturbative expansion of the Hairy 
black hole solution at low orders in perturbation theory. See section |2] for explanation of 
the notation etc. 



B.l Near Field Expansion 

f in (v) _ {y 2 - 1) foV ~ 32) 
j(o,o)\y>- 3e 2 y 4 

fin , v _ {y 2 - 1) (3 (y 4 + y 2 ) e 4 - 96e 2 + 1024) 
J(o,2){y> ~ ^4 

fin _ 32 (y 2 - 1) [27e 4 + 1536e 2 + 384 (3e 2 - 32) log [(l - g) R 2 ] - 22528] ^ 
J(o,4)W- 27e 6 y 4 



f$,o)(v) 



8 {y 2 - 1) (3e 2 (7y 2 - 4) - 64) 
63eV 
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2„,4 



3e z y 

(y 2 - 1) (3eV - 32) 

3y 4 (3 (y 4 + y 2 ) e 4 - 96e 2 + 1024) 

(y 2 - 1) (32 - 3e 2 y 2 ) 2 
2048y 4 (9 (5y 2 - 6) e 4 + 96 (l - lly 2 ) e 2 + 6656) 

3e 2 (y 2 - 1) (3e 2 y 2 - 32) 3 
9e 2 y 6 ^3 (y 3 + y) 2 e 4 - 96 (2y 2 + 3) e 2 + 2048y 2 ) 



+ 



+ 



(y 2 - 1) (3eV - 32) 3 
y 4 (y 8 - 12288 (3e 2 - 32) (3e 2 y 2 - 32)) 

3e 2 (y 2 - 1) (3e 2 y 2 - 32) 3 
32e 2 (40 - 3e 2 ) y 4 



log 



7(y 2 - 1) (32 - 3e 2 y 2 ) 2 



(B.2) 



^2)(y) = 
^(o,4) (v) = 
^(2,0) (v) = 



(6e 2 - 64) 



189e 4 + 4416e 2 + 1536 (3e 2 - 32) log ((l - ^) i? 2 ) - 80896 ' 



18e 5 



2 (3^ 2 + 16) / _ 1_\ 

2i e V W 



y 2 



(B.3) 



& (1,0) 



(v) =1 



C, 2) (y) =« + 



3e 2 



6e 2 y 2 - 128 log (3e 2 - 32) log 



( y 2 -i 

\3e 2 y 2 - 32 



192 log (3e 2 y 2 - 32) 



+ 6 log (3e 2 y 2 - 32) e 2 + 128 log 
+ 64 log 2 (3e 2 y 2 - 32) + 128Li 2 



3e 2 (y 2 - 1) 
3e 2 - 32 

32-3e 2 y 2 
32 - 3e 2 



log (3e 2 y 2 - 32) 



1 



^(3,o) (y) = ^( 150 " 13e ) 

where 

2 (-9e 2 - 192 log (3e 2 - 32) + 288) , 2 (3e 2 - 32 log 2 (32 - 3e 2 ) + 32) 
a = 9? l0g(3) 3? 

+ -^r ~ 18 (e 2 - 32) \og{R) + 61og(e) (3e 2 + 64 log (3e 2 - 32) - 96) 



(B.4) 



(B.5) 
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B.2 Far Field Expansion 



rout i 

Ao,o)l 



-A + r 1 



32 \ 1 
1+ 3? U 



'(0,4)* 
tout 



32e 2 + (1024 - 3e 2 (e 2 + 32)) r 2 



rout j 
J(2,0){ 



3e 4 r 4 
8 (r 4 + 3r 2 + 3) 



9{r 2 + l) c 
1 



[-256 (84r 10 + 463r 8 + 914r 6 + 755r 4 + 193r 2 



189e 2 r 2 (l + r 2 ) 4 
+ 8e 2 (252r 10 + 1261r 8 + 2419r 6 + 2169r 4 + 921r 2 + 99) 



+84 



(3e 2 - 32) r 2 (r 2 + l)* + (32 - e 2 ) r 2 (r 4 + 3r 2 + 3) 



log 



Co)( ? 



1 + r 2 

32 

= 1 + 7^7 



„out 



9(0,4) 
9(2*2) ( 



(r) 



3e 2 y r 2 (l+r 2 ) 2 
9 (r 4 + r 2 + 1) e 4 + 96 (3r 4 + 2r 2 + l) e 2 - 1024 (3 (r 4 + r 2 ) - l) 



8r 2 (r 2 + 3) 
' 9(r 2 + l) 5 



=8 -127e 2 - 



9e 4 r 4 ( r 2 + 1 y 



448 (6r 6 + 20r 4 + 14r 2 + 5) 
(r 2 + l) 4 



+ 2720 



56e 2 (9r 6 + 36r 4 + 70r 2 + 13) 
(r 2 + l) 4 

2 (e 2 (I2r 8 + 57r 6 + 72r 4 + 70r 2 + 13) - 384r 4 (r 4 + 4r 2 + 3)) 

(r 2 + l) 4 



log 



Aout , . 
^(0,0) 



^(0,4) V ) 

A Out I 

(2,2) v 



4 

2 (e 2 (3 



4r 



)-32) 



189e 4 + 4416e 2 + 1536 (3e 2 - 32) log ((l - ||) i? 2 ) - 80896 64 - 6e 



18e 5 

_ e (r 4 + 3r 2 + 3) 9e 2 - 64 
6(r 2 + l) 3 + 42e 
- 75969) e 4 + 96 (23507 - ^) e 2 - 13746176 



+ 



g3 y»2 



26460e 3 

8 (24r 8 + 60r 6 + 20r 4 - 51r 2 - 29) 
9er 2 (r 2 + l) 4 



+ 



e (72r 8 + 228r 6 + 228r 4 + 55r 2 - 35) 



2 (-32 (r 2 + 2) r 4 + e 2 (3r 4 + 8r 2 + 6) r 2 + 64) 
+ — a -log 

3e(r 2 + l) 3 



36r 2 (r 2 + 1) 

„2 



r 2 + 1 



(B.8) 



"(r 2 + l) 2 

2 (-3e 2 + 6 (e 2 - 32) (r 2 + l) log(r) - 3 (e 2 - 32) (r 2 + l) log (r 2 + l) - 32) 
' 3e 2 (r 2 + l) 3 

64r 6 + 260r 4 + 360r 2 - e 2 (9r 6 + 30r 4 + 34r 2 + 13) + 150 
' 63 (r 2 + l) 6 



(B.9) 



C. The soliton at high orders in perturbation theory 
C.l Explicit Results to C(e 17 ) 

As we have mentioned in section ||, the perturbation theory that generates the soliton 
solution as a function of e is straightforward and hence may be automated on Mathematica. 
We have implemented this automation and used it to generate the ground state soliton 
solution to C(e 17 ). For what its worth, we present the resultant explicit formulae for all 
thermodynamical quantities: the mass, the charge and the chemical potential to 0(e). 
Later in this appendix we will also speculate that our solution develops a singularity at the 
origin at a finite value of e. To aid this discussion we also present formulas for f(r = 0) 
and <p{r = 0) to the same order in e. 
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m = 

0.888889e 2 + (1.9737 - 0.170496e 2 ) e 4 + (10.7168 - 1.77184e 2 + 0.0725209e 4 ) e 6 

+ (76.4861 - 18.5347e 2 + 1.48588e 4 - 0.0394005e 6 ) e 8 

+ (624.015 - 198.755e 2 + 23.5941e 4 - 1.23705e 6 + 0.0241682e 8 ) e 10 

+ (5511.63 - 2173.08e 2 + 340.947e 4 - 26.6063e 6 + 1.03262e 8 - 0.0159449e 10 ) e 12 

+ (51307.1 - 24103.6e 2 + 4697.17e 4 - 485.985e 6 + 28.1541e 8 

-0.865871e 10 + 0.0110442e 12 ) e 14 + (495774. - 270273.e 2 + 62898.2e 4 - 8099.8e' 

+623.327e 8 - 28.6648e 10 + 0.729356e 12 - 0.0079209e 14 ) e 16 + C(e 18 ) 



Q = 

0.166667e e 2 + (0.401814e - 0.0364324e 3 ) e 4 + (2.1931e - 0.373206e 3 + 0.0158055e 5 ) e 6 

+ (I5.6491e - 3.87045e 3 + 0.317549e 5 - 0.00864522e 7 ) e 8 

+ (I27.56e - 41.2748e 3 + 4.98583e 5 - 0.266507e 7 + 0.00531978e 9 ) e 10 

+ (I125.66e - 449.55e 3 + 71.5228e 5 - 5.66677e 7 + 0.223607e 9 - 0.00351592e n ) e 12 

+ (I0470.4e - 4972.41e 3 + 980.312e 5 - 102.7e 7 + 6.03018e 9 

-0.188165e n + 0.00243798e 13 ) e 14 

+ (I01107.e - 55636.5e 3 + 13077.2e 5 - 1701.96e 7 + 132.464e 9 - 6.16548c 11 
+ 0.158912e 13 - 0.00174981e 15 ) e 16 + C(e 18 ) 



I 1 = 

- [4. + (-1.52381 + 0.214286e 2 ) e 2 + (-5.87901 + 1.25159e 2 - 0.0652768e 4 ) e 4 

+ (-37.0661 + 10.6372e 2 - 1.00385e 4 + 0.0312297e 6 ) e 6 

+ (-283.701 + 102.563e 2 - 13.7659e 4 + 0.813892e 6 - 0.0179038e 8 ) e 8 

+ (-2410.37 + 1051.44e 2 - 182.027e 4 + 15.6419e 6 - 0.667552e 8 + 0.0113255e 10 ) e 10 

+ (-21860.7 + 11170.5e 2 - 2363.08e 4 + 264.989e 6 - 16.6181e 8 

+0.552789e 10 - 0.00762253e 12 ) e 12 

+ (-207326. + 121451.e 2 - 30325. 8e 4 + 4184.8e 6 - 344.747e 8 + 16.9581e 10 
-0.461289e 12 + 0.00535403e 14 ) e 14 

+ (-2.03127 x 10 6 + 1.34191 x 10 6 e 2 - 386028. e 4 + 63167.2e 6 - 6431. 58e 8 + 417.306e 
-16.8524e 12 + 0.387334e 14 - 0.00387981e 16 ) e 16 + C(e 18 )] 
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f(r = 0) = 

1. - 2.66667e 2 + (-9.04046 + 1.06893e 2 ) e 4 + (-55.7996 + 11.1848e 2 - 0.566798e 4 ) e 6 

+ (-424.503 + 118.131e 2 - 10.9991e 4 + 0.34307e 6 ) e 8 

+ (-3599.49 + 1276.91e 2 - 170.023e 4 + 10.0766e 6 - 0.224439e 8 ) e 10 

+ (-32626.9 + 14049.7e 2 - 2419.13e 4 + 208.257e 6 - 8.96698e 8 + 0.154554e 10 ) e 12 

+ (-309433. + 156626.e 2 - 33000.6e 4 + 3705.4e 6 - 233.893e 8 

+7.87144e 10 - 0.110374e 12 ) e 14 

+ (-3.03237 x 10 6 + 1.76343 x 10 6 e 2 - 438933.e 4 + 60625.9e 6 
-5019.01e 8 + 249.084e 10 - 6.86274e 12 + 0.080996e 14 ) e 16 + C(e 18 ) 

4>{r = 0) = 

e + (2.38095 - 0.206349e 2 ) e 3 + (13.5366 - 2.20766e 2 + 0.0892759e 4 ) e 5 

+ (99.3332 - 23.6891e 2 + 1.86986e 4 - 0.0488621e 6 ) e 7 

+ (825.529 - 258.875e 2 + 30.261e 4 - 1.56287e 6 + 0.0300946e 8 ) e 9 

+ (7388.22 - 2870.91e 2 + 443.951e 4 - 34.1507e 6 + 1.30689e 8 - 0.0199063e 10 ) e 11 

+ (69458.8 - 32195.6e 2 + 6190.33e 4 - 631.96e 6 + 36.1289e 8 

-1.09675e 10 + 0.0138126e 12 ) e 13 

+ (676349. - 364168.e 2 + 83703.e 4 - 10646.1e 6 + 809.226e 8 - 36.7616e 10 
+0.924182e 12 - 0.00991936e 14 ) e 15 + (6.76167 x 10 6 - 4.14712 x 10 6 e 2 
+1.10887 x 10 6 e 4 - 168826.e 6 + 16007.9e 8 - 967.993e 10 + 36.4553e 12 
-0.781807e 14 + 0.00731015e 16 ) e 17 + C(e 19 ) 

C.2 Particular case e = 4 

Of course, the formulae of the previous subsection are not immediately illuminating. In 
order to extract some (tentative) physical conclusions from those formulae, we specialise, 
in this section, to a particular value of e, namely e = 4. At this specific value of e we were 
able to coax Mathematica into producing results upto 0(e 30 ). We will not explicitly list 
our results, but use them to generate some plots that may carry qualitative lessons. 

We are principally interested in the following question: does our solitonic solution 
develop a singularity (and so cease to exist) past a particular critical value of e (or charge)? 
It would seem intuitively that this should be the case; no solitonic solution should exist at 
a mass greater than the 'Chandrasekhar limit' for this configuration. 

Of course it is far from clear that perturbation theory can capture any phenomenon - 
particularly one as interesting as singularity formation - at finite values of e. Nonetheless, 
in this subsection we will investigate the clues that we can glean from our perturbative 
analysis. 

In Figs f| and || we present a plot of f(r) at e = 0.4 and at e = 0.7. We also present 
a plot of the scalar field cp(r) at the origin r = as a function of e. Note that f(r) is 
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f(r) 
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Figure 4: f(r) for e = 0.4 and e = 4 

f(r) 

2r 




Figure 5: f(r) for e = 0.7 and e = 4 

everywhere positive at e = 0.4 while it goes negative near the origin at e = 0.7. Note also 
that the scalar field behaves quite smoothly at the origin at e = 0.4 but shows a pronounced 
peak near the origin at e = 0.7. 

0(r=O) 




Figure 6: <f>(r = 0) as a function of e for e = 4 

We take these results to indicate that the actual solution develops a singularity at 
some value of e between 0.4 and 0.7. Let us use the vanishing of f(r) at the horizon as an 
estimator of the onset of this singularity. In Fig. |7| we plot f(r = 0) as a function of e. 
Note that this graph goes through the origin at e ~ 0.65. 

In Figs |8| and |9| we plot the mass and charge of the solution as a function of e. Note 
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f(r=0) 




Figure 7: f(r — 0) as a function of e for e = 4 

m 

0.35 r 




Figure 9: Charge of soliton as function of e for e = 4 

that these graphs do not show a pronounced peak near e = 0.65. We take this result to 
indicate that the mass and charge of the soliton are finite as we approach the singularity. 
All of the 'conclusions' of this subsection are at best suggestive. A more serious analysis 
of the perturbative expansions described in this appendix could plausibly yield more solid 
indications as to the existence (or otherwise) of a singularity in the solution (rather than 
simply in the perturbative expansion) at e = 0.65. Numerical solutions to the differential 
equations would also likely yield valuable insights here. We leave all such discussions to 
future work. 
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D. Excited State Solitons 



In this section, we will present explicit solution for the solitons obtained by populating the 
first and the second excited state. 

D.l The first excited state soliton 

The solution obtained by populating the first excited state takes the form 



fir) = (r 2 + 1) - 



(r 8 + 5r 6 + 10r 4 + 6) e 2 



2 (r 2 + 1 



,5 



+ -r (e 2 (939978r 20 + 10428693r 18 + 52677075r 16 + 165426459r 14 

34151040 (r 2 + l) 11 v v 

+334808793r 12 + 425064222r 10 + 328198398r 8 + 171552590r 6 + 73128264r 4 

+33080568r 2 + 13962524) - 12 (2568757r 20 + 28256327r 18 + 141281635r 16 

+434422857r 14 + 871830234r 12 + 1127373654r 10 + 916291332r 8 + 520240710r 6 

+186067332r 4 + 111081024r 2 + 26865882)) + O (e 5 ) 



4 



9( r ) = -5TTT + 



1 r 2 (r 6 + 5r 4 - 2r 2 + 12) 



2 



r 2 + l 2(r 2 + l) 7 



e 



4 



(r 2 (e 2 



313326r 18 + 3476231r 16 + 17559025r 14 



11383680 (r 2 + l) 13 
+44264825r 12 + 61767915r 10 + 58915626r 8 + 58567586r 6 + 66442530r 4 
+35744940r 2 + 14941080) - 4 (2568757r 18 + 28256327r 16 + 1419931 15r 14 
+378069945r 12 + 580690770r 10 + 599020422r 8 + 676146702r 6 + 535404870r 4 



+429618420r 2 + 121663080))) + O 




e (r 8 + 5r 6 + 10r 4 + r 2 + 5) 109e 
16(r 2 + l) 5 + 1232 

1 



(e 3 (2266110r 16 - 3679830r 14 - 57389178r 



12 



273208320r 2 (r 2 + l) 11 
-149388558r 10 - 172806150r 8 - 103232855r 6 - 29637025r 4 - 2330873r 2 - 

1065873) + 3e (355291e 2 - 11254468) (r 2 + l) 11 + 12e (-4407480r 16 + 6142752r 14 

+116586624r 12 + 322636512r 10 + 364337655r 8 + 252086890r 6 + 50963352r 4 

+7639614r 2 + 2813617))) + O (e 5 ) 

(2 - 3r 2 ) e e 3 

Mr) = ± + 5 ( e 2 (l308r 12 + 6684r 10 + 13380r 8 + 12637r 6 

y> 2(r 2 + l) 3 9856(r 2 + l) 9V 1 

+5460r 4 + 117r 2 - 710) - 4 (7632r 12 + 41946r 10 + 90252r 8 + 87853r 6 

+41412r 4 - 1053r 2 - 5930)) + O (e 4 ) 

(D.l) 

Here CI is a constant that will be determined by the regularity and normalisability condi- 
tions on <p at one higher order. 
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D.2 The second excited state soliton 

The solution obtained by populating the second excited state is 

16 (r 12 + 7r 10 + 21r 8 - llr 6 + 82r 4 - 18r 2 + 9) e 2 



/(r) = (r 2 + l)- 



45 (r 2 + l) 7 



e 4 



+ ~ (lie 2 (I45257733r 28 + 2191163280r 26 

182614682250 (r 2 + l) 15 v v 

+15436521240r 24 + 71111047280r 22 + 204078916860r 20 + 363116451984r 18 

+447815797000r 16 + 559745609280r 14 + 629383835652r 12 + 760236964032r 10 

+508886827176r 8 + 129434743200r 6 + 41946040800r 4 - 4941526368r 2 

+3675474942) - 128 (809908894r 28 + 12148633410r 26 + 85040433870r 24 

+383937720530r 22 + 1113624190770r 20 + 2086530679662r 18 + 2767907054770r 16 

+3367219159590r 14 + 3524605653426r 12 + 4551299991966r 10 + 1947299375658r 8 

+1248211079850r 6 + 39350124900r 4 + 27778828956r 2 + 14011941561)) + O (e 5 ) 



1 16r 2 (r 10 + 7r 8 - 9r 6 + 85r 4 - 50r 2 + 30) e 



2 



9{r) r 2 + l + 45(r 2 + l) 9 

2 4 

+ — TTf (256 (404954447r 26 + 6074316705r 24 + 42610397025r 22 

182614682250 (r 2 + l) 17 v v 

+146943860245r 20 + 33241 1328795r 18 + 736429945251r 16 + 1692683167175r 14 

+2288373934275r 12 + 2999779439085r 10 + 693256048485r 8 + 761470058349r 6 

+46619547975r 4 + 124153839810r 2 + 57083996970) - lie 2 (l45257733r 26 

+2191163280r 24 + 15436521240r 22 + 48745664240r 20 + 98672175420r 18 

+226204362384r 16 + 539605666600r 14 + 922637944800r 12 + 647631553140r 10 

+673336503840r 8 - 167785201584r 6 + 151437686400r 4 - 3202158960r 2 

+21249708480)) + (e 5 ) 

8 2 /4741e 2 - 228352 e e (8 (5r 4 - 4r 2 + 4) r 4 + l) \ 
A(r) - - + e nn(1 nn„ ^77Z2 + 



90090e 30r 2 30r 2 (r 2 + l) 7 I 

+ e 4 ( CI + (lie 3 (776575800r 24 - 4643959320r 22 

V 973944972000r 2 (r 2 + l) 15 v v 

-49995986040r 20 - 153300661800r 18 - 226524451725r 16 - 196951063595r 14 

-81174009917r 12 - 35347997139r 10 - 47175005605r 8 - 22830663835r 6 

-9008399805r 4 - 109267667r 2 - 82336064) 

+64e (I4151511e 2 - 890269174) (r 2 + l) 15 - 128e (321441 1200r 24 
-18271242990r 22 - 210141141210r 20 - 667533287850r 18 - 983744254350r 16 
-972471806735r 14 

-123178436381r 12 - 380041323207r 10 - 583046341 15r 8 - 160254872155r 6 
-31564145265r 4 - 979195571r 2 - 445134587))) + O (e 5 ) 

(D.2) 
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Here again CI is determined at one higher order. Finally the scalar field in this case is 
given by 



(2r 4 -4r 2 + l)e 2e 3 , 9/ 1S 1fi , A 

(j) r = ^ ^ ~ lie 2 4310r 18 + 27530r 16 + 69735r 14 

(r 2 + l) 4 675675 (r 2 + 1) 12 V V 

+80920r 12 + 29848r 10 - 54222r 8 - 681 10r 6 - 26098r 4 - 6504r 2 + 1149) 

-2 (I141760r 18 + 7624220r 16 + 20257200r 14 + 24337600r 12 + 10397296r 10 

-17335794r 8 - 19274920r 6 - 9051766r 4 - 1800468r 2 + 416883)) + O (e 4 ) 

(D.3) 

E. The First Excited Hairy Black Hole 

In this Appendix we present the results of our construction of the first excited hairy black 
hole metric. We have written a Mathematica programme that allows us to generate the 
corresponding solution for the n th excited hairy black hole at any given value of n. We use 
the same conventions as in appendix^. 

E.l Near Field Expansion 

(y 2 - 1) (eV - 24) 



/(o,2)(y) 

/ ( o,4)(y) 
fRo)(y) 



(y 2 - 1) (e 4 (y 4 + y 2 ) + 304e 2 + 1536) 

8 (y 2 - l) (9Cie 5 + 361e 6 y 2 + 57e 4 (84y 2 + 19) + 576e 2 (26y 2 + 19) + 27648) 
' 9e 6 y 4 

8 (y 2 - 1) (e 2 (l54y 2 - 23) - 212) 
" 231e 2 y 4 

(E.l) 



e 2 y 4 



4 A) {V) = p2( ,_ u y L l2 _ 9d ,s {^eW - We* + eV + if - ^ (2y 4 + 



(y 2 - 1) (e 2 y 2 - 24) 

y 4 (e 4 {y A + y 2 ) - 96e 2 + 3456) 
1 (y 2 - 1) (e 2 y 2 - 24) 2 

t 

e 2 {y 2 - 1) (e 2 y 2 - 24)' 

4-6912e 4 (y 4 + l) + 124416e 2 (y 2 - 4) + 8957952 

8e 2 (712 - 23e 2 ) y 4 
"231 (y 2 - 1) (e 2 y 2 - 24) 2 

(E.2) 
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6 (y 2 - 1) 



Am 12(e 2 -36) (y 2 -l) 



(E.3) 



231ey 2 



'(1,2) 



(y) =^ ( - 288Li 2 (^^) - 12 (e 2 - 72) log(i?) + e 2 (6 y 2 - 5) 



i in 



,(v) 



- 6 (^48 log j - 24 log (e 2 y 2 - 24) + e 2 - 72 j log (e 2 y 2 - 24) 

- 144 log 2 (24^2) + 12 ( e2 " 72 ) log ( e ) " 48 ^ + 456 ) 
5 (71e 2 - 2372) 



(E.4) 

E.2 Far Field Solutions 

f°o%)(r) =r 2 + l 
f out e 2 + 24 



r 



^2 y» 2 



w ^ 1 / 9 / 96 (e 2 - 36) \ ' 



2 (r 8 + 5r 6 + 10r 4 + 6) 



^0)W- 9(r 2 + 1)0 

f(2 U U r ) = " r( ~ 2772r 2 (r 2 + l) (2 (3e 2 - 88) r 12 + 12 (3e 2 - 88) r 10 

{2,2) 1386e 2 r 2 (r 2 + l) 6 V ; V ' 

+ (89e 2 - 2568) r 8 + 5 (23e 2 - 632) r 6 + 80 (e 2 - 24) r 4 + 12 (3e 2 - 88) r 2 

+256) (log(l - ir) + log(l + ir) - 21og(r)) + e 2 (l6632r 14 + 110675r 12 

+309542r 10 + 467241r 8 + 416782r 6 + 179030r 4 + 35422r 2 + 2952) 

- 24 (20328r 14 + 138869r 12 + 394622r 10 + 595703r 8 + 518662r 6 + 240194r 4 



+22250r 2 - 424 



(E.5) 
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r 2 + 1 
e 2 + 24 

g 2 r 2 (j2 _|_ 

e 4 (r 4 + r 2 + 1) + 24e 2 (4r 4 + 3r 2 + l) - 576 (6r 4 + 6r 2 - l) 
" e 4 r 4 ( r 2 + X )3 

2r 2 (r 6 + 5r 4 - 2r 2 + 12) 



9(r 2 + l) 7 



1 



1386e 2 (r 3 + r) 



77 



(r 2 + l) c 



(36 (e 2 - 72) (r 8 + 6r 6 + 3r 4 + 10r 2 



+12) r 4 (log(l - ir) + log(l + ir) - 21og(r)) - e 2 (32r 10 + 182r 8 + 856r 6 
+393r 4 + 210r 2 + 19) + 24 (l08r 10 + 586r 8 + 1292r 6 + 1471r 4 + 222r 2 + 69)) 



- 2567e 2 + 161688 J 



(E.6) 



AOUt I \ 



l (0,0) 

a out ( 
A (0,2)l 

Aout ( \ 
j^out 



r 



e 2 (12 - ^) - 432 



--Ci 



12 e 2 



36) 



e 3 r 2 



AOUt I \ 
A {2,2) V ) 



1 



77e (r 8 + 5r 6 + 10r 4 + r 2 + 5) 



2772 
=C 2 + 



+ 109e 



2544 



(r 2 + l) 5 ^ y 

— (e 2 ((3 (5681r 8 + 30852r 6 + 67813r 4 + 77507r 2 



(E.7) 



11088er 2 (r 2 + 1) 
+39246) r 2 + 29647) r 2 + 1104) + 2772r 2 (r 2 + l) ((5e 2 - 104) r 10 
+8 (3e 2 - 56) r 8 + 5 (9e 2 - 136) r 6 + 40 (e 2 - 8) r 4 
+8 (3e 2 - 56) r 2 + 256) (21og(r) - log (r 2 + l)) 
- 24 (I9673r 12 + 103716r 10 + 217325r 8 + 230143r 6 

424) ) 



+123462r 4 + 10777r 2 
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3r 2 — 2 

(r) =— , 

V ( 1 ' ) V ; 3(r 2 + l) 3 

<*(S)« = 2e2 (r ] +1) 4 ( 3 ( 72 " 5e2 ) ^ + 6 ( e2 " 72 ) ( 3 ^ 4 + r2 ~ 2 ) MO 

-3 (e 2 - 72) (3r 4 + r 2 - 2) log (r 2 + l) - 5e 2 + 456 ) 

1 / (E.8) 

0-* ( r ) = i U ( r 2 + 1) ( 1308r 10 

{3,0)V ' 33264 (r 2 + l) 9 V v 1 v 

+5376r 8 + 8004r 6 + 4633r 4 + 827r 2 - 710) 
+ 4212r 2 - 4 (7632r 8 + 41946r 6 + 90252r 4 

+87853r 2 + 41412) r 4 + 2372o) 

In the above formulae C\ and C2 are constants which will only be determined by 
regularity and normalisability of the scalar field at one higher order. 

F. Thermodynamics in the Canonical Ensemble 

In the following subsections, we will describe the canonical phase diagrams that result 
from a competition between small RNAdS black holes, small hairy black holes and small 
solitons. Everywhere in this section we completely ignore large black holes and large hairy 
black holes. 

In the microcanonical ensemble this was logically justified; large black holes never 
have small mass and charge. However large black holes can (and do) have temperatures 
(and or chemical potentials) comparable to their small counterparts. Consequently the 
phase diagrams we will draw in this appendix do not, in general, represent the true thermo 
dynamical equilibrium of our system at finite temperature and chemical potential. The 
phase diagrams of this appendix should be regarded as formal; their purpose is to help 
us better understand the formal interrelationship between the phases constructed in this 
paper, ignoring all other phases that might exist in the system. 

In this section we will study the interelationship between small black hole, soliton and 
excited black hole phases at fixed charge and temperature. We find it convenient to work 
with the rescaled inverse temperature variable 

^ = St' 

In this section we will assume that /3 and q are small. We also assume that q and f3 2 are 
of the same order, and present all formulae only to leading order in q and /3 2 . 

F.l RNAdS Black Hole 

At any fixed charge, the temperature of a small RNAdS black hole is given, as a function 

W 

^ = 2^(1 - \n 2 ) (F.l) 



of its chemical potential [i = Si by 



-54- 



Small RNAdS BH - determining the 
Domain of existence in the (yS,q) plane 

J* r 



p 




Figure 10: Canonical Ensemble : Small RNAdS BH - For any given charge and temperature 
there are two possible small black holes with two different chemical potentials (Large /i branch has 
a lower free-energy). And for a given temperature, there is a maximum possible charge which is 
attained by the black hole with /i = fi m = •v/3/10. 



In Fig. 10, we present a plot of versus Note that 



2 3 
a < - 
P ~ 2 



(the constraint follows directly from the requirement of positivity of the temperature in 
(j\jD ). Note also that 

"24" 



P - 5 



n l/4 



(F.2) 



This inequality is saturated at at ^ = [i m = J jfc . 

Note that there exist two small black holes (with different values of fi) for any given j3 
that obeys QF.2D The free energy is given as a function on q and fi by 



/ 3 
F sbh = q - (10» + - 



(F.3) 



In Fig. 11 we present a plot of versus /i. Note that 




and the minimum value occurs at /jl = fi m = jq (this is the same value of [i at which 
the temperature curve has a maximum). As is visually apparent from these graphs, the 
RNAdS black hole with the larger value of \i has lower free energy at any fixed (3 and (/(see 
fig. 12). In this appendix we will refer to this solution as the small RNAdS black hole. We 
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Small RNAdS BH - Helmholtz Free energy 
per unit charge vs // 
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Figure 11: Canonical ensemble : small RNAdS BH - Free energy of the small BHs are positive 
and for the large y, branch, the free-energy varies over a bounded domain. 

Small RNAdS BH - Helmholtz Free energy 
per unit charge vs 4n T^Tq 
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Small //^^ 


nch 
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Figure 12: Canonical ensemble : small RNAdS BH - The large n branch always has a lower 
free-energy. 



will completely ignore the free energetically subdominant RNAdS black hole in the rest of 
this Appendix. 



Let us briefly summarise. RNAdS black holes exist whenever the inequality (F.2) is 
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obeyed. Their free energy is given as a function of /3 and q by ( |F.3| ) and (F.l), where we are 
instructed always to choose the larger of the two roots when inverting (F.l). The chemical 
potential of these black holes obey 




F.2 Soliton 



Solitonic solutions exist at all values of f3 and q. At leading order the Free energy and 
chemical potential of the soliton are given by 



Fsoi _ 2ir 
q e 



(F.4) 



Note that ^ < whenever 



32 9 9 9 



e2 ^y x 5 = e ' x 5-- 

It follows that the soliton free energetically dominates that RNAdS black hole whenever 
this inequality is obeyed. At smaller values of e, on the other hand, the RNAdS black hole 
free energetically dominates the soliton at large enough temperatures (but temperatures 
that are small enough to be allowed by ( |F.2j ), i.e. whenever 

3 B 3 

< — < 



2^(3 - 2/4) ~ y/q- 2^(3 - 2tf) 
4 




(F.5) 



We will return to a more detailed comparison of phases below. 



F.3 Hairy Black hole 



Hairy black holes exist if and only if 



In this regime their chemical potential, mass and free energy are given by 



,'4(3e 2 -32) 3 \ _ 16 > 
m = ' 27e6 — r + 3e^ + ' m ' m<L 9 



4 /8(32-3e 2 ) 2 (e 2 -32)\ 2 /9 64^ , ., 



Wq 4(3e 2 - 32) 3 /3 2 \ 16 (32 - 3e 2 ) 4 (21e 4 - 384e 2 + 5120) 4 
~3e + 81e6 J + 1701e 12 ^ 



32 (32 - 3e 2 ) 2 (e 2 - 32) 2 2 (9e 2 - 64) 

■ " ^ -i- *^ y >y ) 

(F.7) 



+ V 63e 7 ^ + 2l7 V + ° ^ V ' ^ 



(where we have listed perturbative corrections, but will only use leading order results in 
what follows). Note that the difference between the free energy of a hairy black hole and 
the soliton is a positive number times /3 2 , so that hairy black holes are free energetically 
always subdominant compared to the soliton. 

It is also interesting to perform a comparison between RNAdS and hairy black holes. 



First notice that the existence ranges (F.6) and ( F.2 ) overlap only when 



9 32 q 

e < — x 5 = ei 
~ 3 2 

Moreover hairy black holes only exist for e 2 > ^ = e 2 . Within this range of e hairy and 
RNAdS black holes in an overlapping region in the /3, q plane. It turns out that RNAdS 
black holes always dominate over hairy black holes in these overlapping regions. 

F.4 Plots of Phase Existence and Dominance 

As we have explained above, the phase diagram of our system depends qualitatively on the 
value of e 2 . In particular, there are three special values of e 2 



3' 1 c 5' 

Recall that RNAdS black holes are always stable - and no hairy black holes exist - for 
e 2 < e 2 . In this regime the only phases of the system are the RNAdS black hole and the 
soliton. It turns out that the RNAdS black hole is free energetically dominant whenever it 



exists. The phase diagram is depicted in Fig. 13 below. 



RNAdS black holes exist only above the line drawn in Fig 13, and give the dominant 
phase when they exist. The soliton dominates the phase diagram elsewhere. 

In the regime e 2 < e 2 < e 2 hairy black holes exist as a phase (below the topmost line 



in Fig. 14) but are never dominant. In this regime the RNAdS black holes exist above the 
bottom most line but free energetically dominate the soliton only in the region between 
the bottom most line and the intermediate line. The soliton is free energetically dominant 
elsewhere. The RNAdS black hole is free energetically dominant over the hairy black hole 
over their common region of existence. 
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1 

Canonical ensemble - Phase diagram - R— vs q 

4?rT 

Case 1 : e 1 < e r 2 

Existence : Small BH in Region A. Soliton everywhere. 

Dominance : Small BH in Region A. Soliton in Region B. 
1 

fi= 

4ttT 




Figure 13: Canonical ensemble : Case 1 : e 2 < e 2 c . No Hairy black hole. Soliton and small 
RNAdS dominate two different regions. 

In the regime e\ < e 2 < e 2 , the phase diagram (shown in Fig.O) is very similar to that 



in Fig. 14, except that there is no intermediate region; the soliton is the thermodynamically 
dominant solution everywhere. Note that in this regime the RNAdS and hairy black hole 
solutions continue to exist as phases; the RNAdS black hole is free energetically dominant 
over the hairy black hole over this region of overlap (i.e., region B in |l5|). 

We turn, finally to the range e 2 > e 2 .. Note that in this range | < = [i m - Recall 
that the chemical potential of a hairy black hole is equal to - at leading order. It follows 
that the RNAdS black hole component of a hairy black hole, in this regime, has /x < fj, m . 
In other words, in this regime, the RNAdS black hole that lies in the centre of a hairy black 
hole is of negative specific heat. In this regime, also, there is no overlap in the existence 
regimes of RNAdS black holes ( |F.2D and hairy black holes ( |F,6| ). At any given value of 
j3 and q we have at most two phases (soliton and black hole or soliton and hairy black 
hole) and the soliton is always the free energetically dominant phase. The phase diagram 
is displayed in Fig. [l6| 

G. Thermodynamics in the Grand Canonical Ensemble 

In the following subsections, we will describe the grand canonical phase diagrams that 
result from a competition between small RNAdS black holes, small hairy black holes and 
small solitons. Everywhere in this section we completely ignore large black holes and 
large hairy black holes. Like our discussion on canonical ensemble given in the previous 
appendix, this is not entirely justified if we are interested in the actual phase diagram 
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1 

Canonical ensemble - Phase diagram - B— vs q 

4ttT 

Case 2 : e t 2 < e 2 < e 2 with e 2 = 3/2 e 2 

Existence : Small BH in Regions A,B & C. 

Hairy BH in Regions B,C & D. Solilon everywhere. 

Dominance : Small BH in Region C. Solkon elsewhere. 
I 

0= 

4nT 




Figure 14: Canonical ensemble : Case 2 : e 2 < e 2 < e\. Hairy black holes appear and over 
regions B and C overlap with the domain of existence of the plain black holes. But, soliton is 
dominant over most of the diagram except for the region marked C where small BHS are preferred. 



of the system. However, the formal diagrams which we present in this appendix are still 
helpful in contrasting the various phases constructed in this paper. 

In this section we will study the thermodynamics of our system as a function of 

R 1 A 4 
AttT e 

G.l RNAdS Black Holes 

As we have seen above, small RNAdS black holes exist only for 




Moreover these black holes also satisfy the inequality 



H " 32(1 - 

This last inequality is automatically obeyed for parametrically small values of f3, of prime 
interest to us in this paper, and so will play no important role in the analysis below. 
Whenever these inequalities are satisfied, we have a unique small black hole. 
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Canonical ensemble - Phase diagram - 8- vs q 

4ttT 

Case 3 : ei 2 < e 2 < c j 2 2 with e 2 = 2 e 2 
Existence : Small BH in Regions A & B. 
Hairy BH in Regions B & C. Solilon everywhere. 
Dominance : Soliton everywhere. 



4jtT 




Figure 15: Canonical ensemble : Case 3 : e\ < e 2 < e?>. Soliton is preferred everywhere over the 
RNAdS or the hairy black hole. 



The various important thermodynamical quantities of small RNAdS black holes, in 
the grand canonical ensemble, are given by 



m 



4(32-3e 2 ) 2 (3e 2 + 32) /3 2 
: 27^ 

64 (e 2 + 32) (3e 2 - 32) 



+ 



9e 5 



16(3e 2 -32) 3 (5e 2 + 32) 
27^ 

/3 2 5/x + 0(/3 5 ,/3 3 ^,<5/i 2 /3) 



P 4 



16 (32 -3e 2 ) 2 /? 2 
9e5 



+ 



256 (3e 2 - 32) 1 



27e 7 



+ 



4 (9 (e 2 - 64) e 2 + 5120) 
9^ 



P 2 5n + ((3 5 ,/3 3 5^,5^/3) 



(G.l) 



4 (3e 2 - 32) V + fl6(32-3e 2 ) 4 \ ^ _ / 6 4 ( 32 3e 2 f\ 



81e 6 

+ (/3 5 ,/3 3 5fi,5fi 2 (3) 



81e 8 



27e 5 



(G.2) 



G.2 Soliton 

The soliton exists for all temperatures but for \i > -. Its thermodynamical quantities are 
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Canonical ensemble - Phase diagram - B— vs q 

4ttT 

Case 4 : e 2 2 < e 2 with e 2 = 1 00 e 2 
Existence : Small BH in Region A . 
Hairy BH in Region C. Soliton everywhere. 
Dominance : Soliton everywhere. 



1 

0= 

4nT 




Figure 16: Canonical ensemble : Case 4 : e| < e 2 . The region of overlap between hairy black 
holes and RNAdS black holes disappears. Soliton is still preferred everywhere. 



given by 



112e 



m 



3 (9e 2 - 64) 



7e 2 



9e 2 
3vr " 



64 



5fi + 



e 2 (2364219e 4 - 47285088e 2 + 244052992) 

1485 (9e 2 - 64) 3 
(I802889e 4 - 39301728e 2 + 215667712) 



6fi 2 + O (<^ 3 ) 



G(jm,T) =— 



Ue 2 



192 - 27e 2 
+ O (<5 M 4 ) 



8/j 2 + O (Sfi 3 ) 

7920 (9e 2 - 64) 3 V ' 

(-1802889e 7 + 39301728e 5 - 215667712e 3 ) \ 
17820 (9e 2 - 64) 3 j 



(G.3) 



It is easy to verify that the solitonic solution always has lower grand free energy (it is 
negative at leading order) than the RNAdS black hole (the free energy is positive at leading 
order ) within the validity of perturbation theory. Consequently the system undergoes a 
first order phase transition from the RNAdS black hole to the solitonic phase as \i is raised 
above 4 



G.3 Hairy Black hole 

Hairy black holes exist whenever 



bjx_ 8 (3e 2 - 
]P ~ 9^7 



32) : 



(G.4) 
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Their thermodynamical quantities are given by 

4 (252e 7 5/x + (32 - 3e 2 ) 2 (27e 4 - 576e 2 + 5120) /3 2 ' 

m =- 27^-64) ' + W"> ** 

_21e^-8(32-3e 2 ) 2 (e 2 -32)^ + o(/34;/32 ^^ ) 



3e 5 (9e 2 - 64) 

4 (3e 2 - 32) 3 /? 2 16 (32 - 3e 2 ) 4 (27e 6 - 696e 4 + 10752e 2 - 57344) /3 4 
81e6 + 243e 12 (9e 2 - 64) 



2 (21e 7 fy 2 - 16 (32 - 3e 2 ) 2 (e 2 - 32) f3 2 5n 
9e 5 (9e 2 - 64) 



(G.5) 



It is easily verified that (within perturbation theory) hairy black holes are free en- 
ergetically subdominant compared to solitons in their common domain of existence. On 
the other hand, it may be checked that they are free energetically dominant compared to 
RNAdS black holes, where the solutions coexist. 

G.4 Phase Diagrams 

The phase diagram of our system is very simple when e 2 < e 2 . Hairy black holes don't 
exist. The two phases that do exist - RNAdS black holes and the soliton - never coexist at 
the same [i and f3. The RNAdS black hole dominates when it exists; the soliton dominates 
when it exists. The phase diagram is sketched in Fig.[H] below. 

The phase diagram is more interesting when e 2 > e 2 . As we have mentioned above, 
the system undergoes a first order phase transition from the black hole to the solitonic 
phase as fi is raised above ~. The hairy black hole phase is always subdominant compared 
to the soliton, but free energetically dominates the black hole when both exist. The phase 



diagram is depicted in Fig. 15 below. Note that the black hole and hairy black holes phases 
are identical, where the hairy black hole is first created. In the absence of the solitonic 
solution, consequently, our system would have undergone a second order transition from 
the RNAdS black hole to the hairy black hole phase upon raising /i. 

H. Notation 
H.l Basic Setup 

Throughout this paper, we work in asymptotically (global) AdSs spacetimes with a bulk 
metric g, a bulk charged scalar field (j) and a bulk gauge field with a Lagrangian 



5 =8^/^ 



\ (K[g\ + 12)- I 



(TZ[g] + 12) - -J» V T» V - \D^\ 2 ~ ™\\<S>? 



(H.l) 



T^v = V^A U - V u Afj, and D^(f> = V M (/> - ieAtf 



where G5 is the Newtons's constant and the radius of AdSs is set to unity. This implies 
that the bulk cosmological constant is taken to be A5 = —6. 
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Figure 17: Grand canonical ensemble : Case 1 : e 2 < e 2 c . No hairy black holes. Solitons and 
RNAdS black holes never overlap and they dominate their respective domains of existence. 



The radial co-ordinate of AdSs is denoted by r with the boundary of AdSs being at 
r = oo. For solutions with the horizon, the outer horizon is taken to be at r = R. The 
Schwarzschild-like temporal co-ordinate is denoted by t. Sometimes, we find it convenient 
to work with rescaled co-ordinates y = r/R and r = t/R, especially in the near-field 
expansion at small radius (r << 1). In Appendix [A], we shift to Eddington-Finkelstein like 
co-ordinates, with an Eddington-Finkelstein(EF) time co-ordinate denoted by v. 

Since throughout this paper we work with spherically symmetric solutions, we will 
leave the co-ordinates parametrising S 3 implicit. Mostly, we work in a gauge where the 
bulk fields take a form 

ds 2 = -f(r)dt 2 + g{r)dr 2 + r 2 dQl 



A t = A(r) 







(H.2) 



In Appendix |A], we work with the EF co-ordinate metric for a small charged black hole 
(see below). 

The charge of the scalar field (j) is denoted by e and its mass by m^. We take = 



for most of the paper except in section . By the standard rules of AdS/CFT, the dual 
boundary operator 0$ has a scaling dimension 



A r 




2 + JA + m 2 



d=4 
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Grand Canonical ensemble - Phase diagram - 8- vs u 

4jtT 

Case 2 : e 2 < e 1 with e 2 = 1.51 e c 2 

Existence : Soliton in Region B,C,D & E. 

Small BH in Regions A,B & C. Hairy BH in Region D & E. 

Dominance : Small BH in Region A. Soliton elsewhere, 
i 

H= — 

4zrT 




P. 



Figure 18: Grand canonical ensemble : Case 2 : e 2 c < e 2 . Small BH is preferred region A and 
soliton is preferred elsewhere. 

This is also the energy of the lowest <j> mode in vacuum A0IS5. For the case = 0, 
Ao = 4. The other spherically symmetric modes of <f) (dual to the descendants d 2n O,p ) 
have an energy 

A n = A + 2n = 2 + y/ITm^ + 2n 

For the case = 0, A n = 4 + In. The covariant derivative acting on <p is denoted by 
Dfj, = V M — ieA^. The symbol D 2 is used to denote the covariant Laplacian. 

H.2 Thermodynamic quantities 

We now turn to notations involving thermodynamic quantities. First, we omit a factor 
of G5 1 from all our extensive quantities in order to simplify our expressions. With this 
understanding, we will denote the ADM mass by M, the charge of a solution by Q and its 
entropy by S. We often find it convenient to work with a rescaled mass m and a rescaled 
charge q which are related to the actual mass M and charge Q via the relations 

Q = -q and M = —m (H.3) 
2 8 

We use F = M—TS to denote Helmholtz free-energy appropriate to the canonical ensemble 
and G = M — TS — fiQ to denote the grand potential appropriate to the grand-canonical 
ensemble. Coming to the intensive variables \x represents the chemical potential and T 
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represents the temperature both of which are defined by the first law 

dM = TdS + udQ 

We sometimes find it convenient to work with the 'rationalised' inverse temperature j3 = 
(47rT) -1 and the chemical potential excess over the super-radiant bound Sfi = /z— A n /e. 

H.3 Double expansions 

Quantities in this paper are often expressed as double expansions about two parameters - 
first is the parameter e which is the amplitude of the leading normalisable mode in cf). Under 
AdS/CFT, roughly e ~ < >, the boundary expectation value of O^, the operator dual 
to <p. The second parameter is the outer horizon radius R of the small charged blackhole 
at the core of the hairy blackhole. Further, our solutions are often expressed in terms of a 
matched asymptotic expansion with a near-field expansion at small radius (r << 1) and a 
far field expansion far away from the horizon (r >> R) matched at their common domain 
of validity. We use the superscripts in and out to denote these two expansions respectively. 
Many formulae in this paper involve the coefficients in this expansion which are defined 
via 



f(r) 

g{r) 

A(r) 



Near field(r « 1) : /*» = £^=0 f£ = EZo ^ ET=o 'ft,2k 
Far field(r » R) : /"* = £~ = E~ o e 2 " £~o 



.fin 


— 2^n=0 6 


j out 




. g in 




gOUt . 


- Z^n=0 fc 


. A in 




AOUt 


v-voo 


. ^in 


v-\oo 


^out 


- V°° f 

— Z^n=0 e 



2n„in v^oo 2n v^oo r>2k n in 

y-ln — l^n=0 t l^k=0 11 y2n,2fc 

ln„out v^oo 2n v^oo r>2k n out 

i/2n — 2_^n=0 e 2^fc=0 il y2n,2fe 

-2n Ain _ \poo 2n v^oo r>2k /Jin 
: ^2n — 2^n=0 fc 2^fc=0 21 ^2n,2fc 
2n Aout _ v^oo 2n v^ 00 E>2fc /tout 
^2n — Z^n=0 fc Z^A:=0 il ^2n,2k 

lield( r « J ! : , - E^o = En= ^ 2n+1 E^o^SU 

2n+,2fc 

(H.4) 



Far field(r » R) : ^ = E^o^W = En=o e2 " +1 Er=o i? 2 *-*^ 



In a similar vein, one can define an expansion of the covariant Laplacian 

= f Near field(r « 1) : (!*»)*> = E~ e 2n (Z? 2 )|™ = £~ ^ E^o ^WSU 
" \ Far field(r » fl) : = E~o ^(^S? = E~=o ^ E^O a 2 *^)^ 

Further, the boundary value of the gauge field or the chemical potential has the double 
expansion 

oo oo oo 

a = lim A(r) = V e 2n ^ 2n (R) = V e 2 " V fl 2 W,2fc 

n=0 n=0 fe=0 
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